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Abstract 

CFD-Based  Approximation  Concepts 
for  Aerodynamic  Design  Optimization, 

With  Application  to  a  2-D  Scrainjet  Vehicle 

by  Peter  D.  Me  Quade 

Chairperson  of  the  Supervisory  Committee:  Professor  D.  Scott  Eberhardt 

Department  of  Aeronautics  and  Astronautics 


This  dissertation  investigates  the  application  of  approximation  concepts  to  aero¬ 
dynamic  optimization.  Such  techniques,  which  are  gaining  popularity  in  structural 
optimization,  offer  the  potential  of  providing  the  accuracy  of  a  high-fidelity  “detailed” 
analysis  model  at  greatly  reduced  computational  cost.  This  is  because  the  detailed 
model  is  used  only  to  “fine-tune”  an  approximate  model  which  is  then  used  in  the 
optimizer.  The  test  problem  treated  is  the  design  optimization  of  a  2-D  scramjet  vehi¬ 
cle  flying  at  Mach  6.0  at  30  km  altitude.  The  objective  function  is  net  thrust.  The 
following  approximation  concepts  are  used;  the  Taylor  series  approximation  to  wall 
pressures  and  inlet  plane  flow  properties;  and  Haftka’s  Global-Local  Approximation 
applied  to  the  same  variables.  The  performance  of  these  techniques  is  compared  to 
that  for  optimization  using  CFD  alone.  Cost  reductions  ore  quantified. 

It  is  shown  that  modifications  must  be  made  to  the  formulation  of  the  approxima¬ 
tion  concepts  as  they  are  used  in  structural  optimization,  due  to  the  changing  grid 
geometries  required  by  the  CFD  solver.  All  correction  factors  for  the  approximation 
concept  are  applied  net  to  the  CFD  grid  points,  but  to  a  constant,  dense,  nondiraen- 
sionaUzed  “correction  point  grid”,  which  does  not  change  as  the  CFD  grid  changes.  It 
is  also  shown  that,  in  areas  where  discontinuous  phenomena  are  not  important  (such 
as  in  the  scramjet  nozzle),  the  approximation  concepts  can  be  successfully  used,  after 
this  modification  is  made.  Optimizations  of  the  nozzle  region  show  that  all  the 


approximation  concepts  result  in  a  68%  reduction  in  the  number  of  calls  to  the  CFD 
solver. 

In  regions  dominated  by  shock  impingements  (such  as  the  forebody/inlet),  it  is 
found  that  approximation  concepts  applied  to  point  properties  cannot  be  used  as  they 
currently  are  in  structural  optimization,  due  to  the  effects  of  shock  movement  during 
correction  factor  calculations,  and  due  to  artifacts  of  the  CFD  solver,  such  as  shock 
smearing.  In  fact,  even  though  the  CFD  and  the  (uncorrected)  approximate  models 
optimize  to  very  nearly  the  same  design,  the  Taylor  series  and  GLA  fail  to  do  so. 
However,  application  of  the  GLA  to  the  integrated  objective  function  (net  thrust)  with 
zeroth-order  correction  factors,  is  successful. 

To  lay  the  groundwork  for  future  investigation,  a  method  of  improving  the  behav¬ 
ior  of  the  point-property  GLA  in  the  presence  of  shock  impingements  is  developed 
and  tested.  This  involves  using  “floating”  pre-  and  post-shock  coordinate  systems  for 
each  wall  surface.  The  result  is  a  dramatic  reduction  in  the  erratic  behavior  of  the 
GLA.  This  technique  may  form  the  basis  of  a  generally-applicable  GLA  technique  for 
aerodynamic  optimization. 
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Chapter  1 


Introduction 


1.1  Background 

The  dramatic  improvements  in  computing  power  and  computational  techniques 
in  recent  years  have  made  computational  fluid  dynamics  (CFD)  a  valuable  tool  for 
analyzing  complicated  internal  and  external  flows.  CFD  has  begun  to  see  use  in  air¬ 
craft  design,  including  the  wing  designs  of  the  next  generation  of  commercial  trans¬ 
port  aircraft,  and  the  preliminary  design  of  the  National  Aerospace  Plane  (NASP). 
Indeed,  for  this  latter  program,  CFD  provides  the  only  means  of  gathering  detailed 
flow  data  for  the  high-speed  portion  of  the  flight  envelope. 

As  CFD  continues  to  mature  and  to  become  more  reliable  and  cheaper  to  use,  it  is 
natural  that  designers  will  seek  to  incorporate  CFD  with  automated  design  synthesis 
(optimization)  methods.  Such  techniques  have  been  used  successfully  in  structural 
design  in  recent  years.  (Ref’s.  7,  11  and  14.)  They  have  also  been  applied  to  aero¬ 
dynamic  design  problems,  but  using  only  relatively  simple  flow  analysis  techniques. 
(Ref’s.  6,  10  and  11.)  The  high  cost  of  using  CFD  as  an  analysis  model  has  gener¬ 
ally  precluded  its  direct  use  in  vehicle  design  opumization,  which  requires  repeated 
analyses  of  the  objective  function  and  constraints. 

Recent  developments  in  structural  optimization  techniques  have  focused  on 
ways  to  reduce  the  cost  of  using  “detailed”  analysis  models  in  optimization.  One 
such  promising  area  of  research  is  “approximation  concepts”.  (See  Ref.  7.)  These 
concepts  are  based  upon  replacing  die  original  problem  with  a  simpler  (and  much 
cheaper),  “approximate”  problem  which  is  used  in  the  actual  optimization  process. 
The  fidelity  of  such  an  approximate  problem  can  be  enhanced  by  using  information 
from  a  limited  number  of  analyses  from  the  detailed  model  to  “fine  tune”  the  approxi¬ 
mate  model  before  it  is  used  by  the  optimizer. 

A  simple,  and  now  commonly-used  example  of  approximation  concepts  involves 
forming  a  first-order  Tiylor  series  approximation  to  the  physical  propenies  of  interest 
(point  stresses,  for  example),  in  terms  of  the  design  variables.  The  coefficients  in  the 

Tiylor  series  come  from  applying  the  detailed  analysis  model  at  a  baseline  design. 
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Naturally,  such  an  approximation  is  only  guaranteed  to  be  valid  in  the  vicinity  of  the 
baseline  design.  For  this  reason,  Taylor  series  approximations  are  referred  to  as  “local 
approximations.”  These  approximations  are  implemented  in  a  design  optimization  by 
performing  a  complete  optimization  within  certain  “move  limits”  about  the  baseline, 
and  then  forming  a  new  Taylor  series  approximation  about  tfie  optimum  design,  estab¬ 
lishing  new  move  limits,  and  repeating  the  process.  The  detailed  model  is  never  actu¬ 
ally  used  in  the  optimizer. 

An  approximation  concept  which  is  potentially  less  limited  than  the  Taylor 
series  approach  is  Haftka’s  “Global-Local  Approximation”  (GLA),  in  which  a  simple, 
but  globally-applicablc  analysis  model  is  “fine-tuned”  by  data  from  a  few  analyses 
from  the  detailed  model.  (Ref.  8.)  The  fine-tuned  approximate  model  is  then  submit¬ 
ted  to  the  optimizer.  The  GLA  is  implemented  by  forming  correction  factors  which 
relate  the  detailed  model’s  results  to  the  approximate  model’s  results  at  a  baseline 
design.  Such  a  correction  factor  may  itself  take  the  form  of  a  first-order  Taylor  series 
approximation  to  the  ratio  of  detailed/approximate  model  results.  Thus,  the  GLA 
would  be  expected  to  provide  very  accurate  results  near  the  baseline  design,  but  retain 
much  better  accuracy  than  the  Taylor  series  approach  as  the  design  leaves  the  vicinity 
of  the  baseline. 

Although  GLA’s  have  been  applied  to  some  structural  problems  (Ref’s.  7  through 
9),  they  have  only  very  mcently  been  applied  to  aerodynamic  optimization.  Ref.  10, 
which  is  perhaps  the  first  published  example,  was  published  in  January,  1992.  GLA’s 
have  not  heretofore  been  applied  using  CPD  as  the  detailed  model. 

The  cost  savings  resulting  from  the  successful  application  of  approximation  con¬ 
cepts  can  be  great  This  is  increasingly  true  as  the  number  of  design  variables 
increases.  For  instance,  a  problem  with  100  design  variables  and  which  is  solved 
using  a  detailed  model  for  all  artalyses  would  require  more  than  100  calls  to  the 
detailed  model.  A  successful  approximation  concept  might  reduce  this  to  10  calls  or 
less  (and  hundreds  of  calls  to  the  much  faster  approximate  analysis  model). 

However,  there  are  potentially  serious  problems  inherent  in  applying  approxima¬ 
tion  concepts,  as  they  have  been  formulated  for  structural  problems,  to  aerodynamic 
optimization.  For  example,  the  presence  of  shock  waves,  transition  to  turbulence, 
and  Bow  separation/reattachraent  are  discontinuous  phenomena  which  could  interfere 
with  the  linear  extrapolations  used  in  the  fine-tuning  process,  particularly  if  these  phe¬ 
nomena  move  as  the  design  is  changed.  In  addition,  artifacts  of  the  CFL)  solution 
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process  (such  as  shock  smearing  and  pre-shock  oscillation)  may  influence  the  hne  tun¬ 
ing  and  introduce  errors  in  calculating  search  directions.  This  is  complicated  by  the 
fact  that  the  detailed  and  approximate  models  may  not  capture  shocks  or  other  phe¬ 
nomena  equally  “crisply”,  which  can  lead  to  erroneous  correction  factors.  Also,  the 
structural  applications  of  approximation  concepts  have  so  far  only  dealt  with  fixed- 
grid  geometries.  For  CFD  to  be  used,  however,  the  grid  must  be  allowed  to  vary  as 
the  geometry  and/or  flow  phenomena  change,  even  to  the  point  of  adding  or  subtract¬ 
ing  grid  points.  This  precludes  a  straightforward  assignment  of  correction  factors  to 
grid  points  as  was  used  in  previous  structural  optimizations. 

There  are  other  problems,  not  related  to  the  approximation  concepts,  but  which 
are  unique  to  applying  nonlinear  optimization  to  aerodynamic  design.  First,  it  is  not 
yet  clear  if  the  sudden  appearance  of  shock  waves,  separation,  or  ether  phenomena  on 
the  aircraft  surface  might  introduce  more  discontinuity  in  such  measures  as  lift,  thrust, 
and  drag  than  can  be  handled  by  a  gradient-based  optimizer.  Also,  as  has  been  shown 
in  structural  optimization,  the  choice  of  design  variables  can  have  a  great  impact  on 
whether  the  optimization  leads  to  the  proper  design.  But  in  aerodynamic  vehicle 
design,  additional  complexity  is  introduced  by  the  fact  that  certain  defining  geome¬ 
tries  of  a  design  are,  in  practice,  not  changed  independently,  but  are  changed  in  con¬ 
cert.  For  example,  a  transonic  swept  wing  might  use  twist  to  maintain  a  desirable  lift 
distiibvtion.  As  the  sweep  is  increased,  the  twist  should  also  be  increased.  One  could 
devise  one  or  more  sets  of  design  variables  which  would  make  this  connection  auto¬ 
matically.  However,  it  is  not  always  clear  at  the  outset  which  set  of  design  variables 
will  work  best  in  an  optimizer. 


1.2  Problem  Statement 

It  is  the  purpose  of  this  research  effort  to  investigate  the  application  of  CFD- 
based  approximation  concepts  to  aerodynamic  design  optimization.  The  test  case 
selected  is  the  design  of  a  2-dimcnsionaI  scramjet-poweied  vehicle  in  an  inviscid,  per¬ 
fect  gas  flow.  (See  Fig.  1.1.)  The  flight  conditions  are:  =  6.0,  h  =  30  km.  Only 

the  lower  surface  of  the  vehicle  is  treated,  since  this  is  where  all  the  dominant  phe¬ 
nomena  occur.  This  test  case  is  illustrative  of  all  the  problems  mentioned  above.  Spe¬ 
cifically,  it  has  legions  dominated  by  relatively  isentropic  flow,  and  regions  dominated 
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by  shock  waves  which  impinge  upon  vehicle  walls.  It  requires  grid  clustering  which 
must  change  as  the  design  is  changed.  And  the  design  is  complicated  enough  that  the 
best  selection  of  design  variables  is  not  obvious  in  advance.  Finally,  the  design  prob  ¬ 
lem  is  significant  and  timely  in  its  own  right,  since  scramjet  vehicle  design  is  cur¬ 
rently  being  vigorously  pursued  in  support  of  the  NASP  program. 

Both  CFD-based  Taylor  series  and  GLA’s  were  formulated,  and  the  unique  prob¬ 
lems  inherent  in  applying  these  approximation  concepts  to  aerodynamic  optimization 
were  investigated.  New  techniques,  not  required  in  structural  optimization,  needed  to 
be  developed  for  this  application.  The  approximation  concepts  were  then  applied  to 
the  optimization  process,  and  their  behavior  and  cost  savings  relative  to  straightfor¬ 
ward  optimizations  using  CFD  alone  were  observed.  The  optimization  package  u.sed 
was  “Design  Optimization  Tools”  (DOT).  See  Ref.  12  and  Appendix  A. 

Because  of  the  inherently  different  natures  of  the  flows  in  the  forebody/inlet 
region  (which  is  dominated  by  oblique  shock  waves)  and  in  the  nozzle/afterbody 
region  (which  is  primarily  an  isentropic  expansion  process),  these  optimization  tech¬ 
niques  were  applied  to  each  of  these  regions  separately.  Thus,  problems  arising  from 
the  various  flow  phenomena  were  observed.  It  was  found  that  such  effects  have  a 
great  impact  on  the  ease  with  which  the  approximation  concepts  can  be  applied. 

In  order  to  provide  a  challenging  test  of  approximation  concepts  in  the  presence 
of  shocks  which  move  as  the  design  changes,  the  forebndy/inlet  optimizations  were 
performed  allowing  arbitrary  shock  configurations  there.  That  is,  the  problem  was  not 
limited  to  one  with  “perfectly-placed”  shocks,  with  the  nose  and  ramp  shocks  imping¬ 
ing  on  the  cowl  lip,  and  the  resulting  cowl  lip  shock  hitting  the  inlet  comer.  This  for¬ 
mulation  of  the  problem  highlighted  the  difficulties  one  would  encounter  in  applying 
GLA's  to  more  general  problems  with  moving  shocks,  such  as  in  transonic  wing 
design. 


1.3  Original  Contributions 

This  research  effort  is  the  first  application  of  tlie  GLA  technique  in  which  CFD  is 
the  detailed  analysis  model.  As  such,  this  research  identifies  general  procedures  for 
using  CFD  with  GLA’s.  It  also  identifies  problems  which  will  be  encountered,  not 
only  with  the  design  problem  at  Jiand,  but  in  all  aerodynamic  design  problems  where 
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dissontinuous  phenomena  (such  as  shock  waves  or  transition  to  turbulence)  exist  and 
move  as  the  design  is  changed.  Indeed,  it  is  shown  that  the  GLA.  as  currently 
employed  in  structural  optimization,  cannot  be  straightforwardly  employed  when  such 
phenomena  are  present.  Proposed  solutions  to  these  problems  are  investigated,  and 
areas  of  further  work  are  recommended. 

This  work  identifies  the  general  problem  that  changing  grids  cause  in  any  approx¬ 
imation  concept.  Although  this  problem  had  been  observed  in  a  different  type  of  opti¬ 
mization  approach  (Ref.  3),  it  was  not  explicitly  solved  there,  but  was  “worked 
around"  by  changing  the  objective  function.  In  this  research  effort  a  direct,  general 
solution  to  this  problem  has  been  found  and  successfully  applied. 

This  is  also  the  first  application  of  GLA*s  to  the  aerodynamic  optimization  of  a 
scramjet  vehicle  design.  It  provides  new  insights  into  the  selection  of  approximate 
models  and  design  variables  for  vehicle  design  optimization. 

Approximation  concepts  hold  great  promise  for  dramatically  reducing  the  cost  of 
aerodynamic  design  optimization  using  CFD.  Research  such  as  this  work  is  an  impor¬ 
tant  part  in  the  effort  to  make  CFD-based  optimization  techniques  a  part  of  the  aircraft 
designer’s  repertoire. 


Fig.  1.1  2-D  scramjet  vehicle  geometry. 


Chapter  2 


Approximation  Concepts 
2.1  Introduction 

The  field  of  approximation  concepts,  as  currently  used  in  structural  optimization, 
has  produced  techniques  which  can  be  used  to  substantially  reduce  the  number  of 
costly  function  evaluations  in  an  optimization  process.  The  basic  principle  of  approx¬ 
imation  concepts  is  to  replace  the  original  problem  with  an  approximate  problem, 
which  is  simpler  and  less  costly  to  solve.  An  approximate  problem  solver  is  then 
implemented  for  function  evaluations.  It  is  this  approximate  model  which  is  then 
submitted  to  the  optimizer:  the  original,  “detailed”  model  is  never  used  directly  by 
the  optimizer.  It  is  possible  to  improve  the  fidelity  of  the  approximate  model  by 
using  information  from  a  limited  number  of  analyses  from  the  detailed  model  to  “fine 
tune”  the  approximate  model.  The  focus  of  this  research  is  on  this  type  of  approxima¬ 
tion  concept. 

For  the  purposes  of  this  research,  the  detailed  model  is  the  CFD  solver.  Several 
approximate  models  were  used  with  the  CFD  solver  to  form  the  bases  of  the  various 
approximation  concepts. 

Approximation  concepts  can  be  applied  to  “point”  properties,  such  as  wall  pres¬ 
sures,  or  to  “integral”  properties,  such  as  thrust,  drag,  or  net  thrust.  If  the  objective 
function  is  itself  an  integral  property,  the  use  of  point  properties  for  the  approxima¬ 
tion  concept  may  capture  nonlinearities  in  the  objective  function  that  the  use  of  inte¬ 
gral  properties  might  not.  Thus,  for  most  cases  in  this  research,  the  approximation 
concepts  were  applied  to  point  properties  (wall  pressures  and  inlet  plane  flow  proper¬ 
ties). 

This  chapter  describes  the  two  approximation  concepts  used  in  this  project:  the 
Taylor  series  “local  approximation”  and  the  Global-Local  Approximation  (GLA).  It 
also  presents  details  of  how  these  concepts  are  implemented  in  the  2-D  scramjet  aero¬ 
dynamic  optimization  problem.  There  are  special  difficulties  encountered  in  applying 
them  to  aerodynamic,  versus  structural,  optimization  problems.  These  difficulties, 
and  insight  into  their  solutions,  are  presented  as  well. 
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2.2  Local  Approximations 

One  often-used  method  of  reducing  the  number  of  costly  calls  to  the  detailed 
analysis  program  in  an  optimization,  while  still  retaining  some  of  the  accuracy  of  that 
detail  is  to  use  “local  approximations**  (using  Haftka’s  nomenclature,  Ref.  8.).  Such 
approximations  extrapolate  the  results  of  a  detailed  analysis  done  at  a  baseline  design 
point,  and  are  therefore  valid  only  in  the  immediate  vicinity  of  that  design  point. 
The  most  common  type  of  local  approximation  is  a  Taylor  series  expansion  of  the  rel¬ 
evant  properties  about  the  latest  baseline  de.sign. 

A  first-order  Taylor  series  approximation  for  a  point  property  can  be  formed  as 
follows.  Let  X  be  the  vector  of  design  variables  (such  as  nose  angle,  compression 
ramp  angle,  etc),  and  let  f(X)  be  the  property  of  interest  (say  the  pressure  at  a  particu¬ 
lar  point  on  a  wall).  Define  X*’  to  be  the  latest  “baseline**  design;  that  is,  the  latest 
design  for  which  we  have  results  from  the  detailed  model.  We  then  have 

f(X)-f^(X®)+Vf,AX,  (2.1) 

where:  f^  is  the  value  of  the  property  found  using  the  detailed  model;  V  f^  is  a  row 
vector  containing  the  elements  of  the  gradient  of  f^  with  respect  to  the  design  vari¬ 
ables;  and  A  X  s  X  -  X" .  Note  that  every  point  at  which  wall  pressures  are  to  be  eval¬ 
uated  has  its  own  Taylor  series  approximation  to  f(X). 

Optimization  using  local  approximations  proceeds  as  follows.  (See  Fig.  2.1.) 
The  detailed  model  is  used  to  find  the  baseline  values  of  f^  and  and  V  fj.  A  complete 
optimization  (or  “iteration*’)  is  then  performed,  following  a  technique  such  as  that  out¬ 
lined  in  the  Appendix.  Since  the  approximation  is  only  valid  in  the  vicinity  of  X",  we 
apply  “move  limits*’,  which  restrict  AX  so  that  no  design  variable  changes  by  more 
than  a  given  percentage  during  the  iteration.  When  the  iteration  is  completed,  a  new 
detailed  analysis  is  performed,  and  new  coefficients  for  the  approximation  are  calcu¬ 
lated.  The  process  continues  until  convergence  is  reached. 

Local  approximations  have  proven  useful  in  structural  optimization  problems 
(Ref’s.  8  and  9),  but  have  only  recently  been  used  in  aerodynamic  optimization  pro- 
lems.  (Ref.  4.)  First-order  Taylor  series  were  employed  in  several  of  the  optimiza¬ 
tions  done  in  this  research,  as  explained  in  Chapters  4-6. 


Detailed  Analysis 


Fig.  2.1  General  application  of  an  approximation 
concept 
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2.3  Global-Locsil  Approximations 

A  method  of  adding  some  degree  of  global  accuracy  to  a  local  approximation  is 
to  use  a  “Global-Local  Approximation”  (GLA),  Sec  Ref’s  8  and  9.  In  this  technique, 
we  use  f«(X) ,  the  result  from  a  (globally-applicable)  approximate  model  of  f(X),  and 
apply  a  linearly-varying  correction  factor  p  (X) ,  which  “corrects”  f^  to  more  closely 
resemble  the  results  of  the  detailed  model,  f^. 

/(X)-p{X)/,(X)  (2.2) 

where 


P(X) 


/,,(X) 

/„(X) 


P(X°)  +  (Vp(X''))AX 


and,  by  chain  rule  differentiation, 


(2.3) 


Vp(X°) 


(2.4) 


Eq.  (2.3)  shows  that  p  (X)  is  represented  by  a  first-order  Taylor  series  about  the 
latest  design  point  for  which  detailed  model  data  is  available.  At  the  conditions  of  the 
latest  design  point,  the  GLA  representation  of  f(X)  is  equal  to  the  detailed  model 
results.  This  is  the  local  nature  of  the  GLA.  Note  that  the  calculation  of  V  p  (X”) 
requires  calculation  of  the  various  elements  of  both  V  (X°) ,  which  can  be  costly,  and 
V  /,  (X®) ,  which  is  much  less  so. 

The  GLA  method  is  applied  in  the  same  general  fashion  as  the  local  approxima¬ 
tion  was.  ( Fig.  2.1.)  For  the  point-property  GLA  applied  to  the  2-D  scramjet  vehicle, 
p’s  are  determined  for  all  surface  pressures,  and  for  the  inlet-plane  flow  properties, 
based  on  the  initial  design.  A  complete  optimization  (or  “iteration”)  is  conducted, 
using  only  the  approximate  model  with  corrections  applied  as  per  Eq.  (2.2).  Since 
the  GLA  is  only  valid  in  the  vicinity  of  X®,  we  apply  “move  limits”,  which  restrict  AX 
so  that  no  design  variable  changes  by  more  than  a  given  percentage  during  the  itera¬ 
tion.  Upon  completion  of  the  iteration,  new  p’s  are  determined,  new  move  limits  are 
established,  and  a  new  iteration  is  conducted.  The  process  continues  until  conver¬ 
gence  is  reached.  GLA’s  were  used  for  several  of  the  optimizations  used  in  this 
research  project,  as  explained  in  Chapters  4-6. 
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2.4  Applying  Approximation  Concepts  to 
Aerodynamic  Optimization 

A  special  problem  is  posed  when  applying  these  approximation  concepts  to 
CFD-based  aerodynamics  problems,  which  necessarily  involve  using  CFD  grids 
which  add  or  delete  grid  points  as  the  design  variables  change  (to  accommodate  a 
lengthened  cowl,  or  a  change  in  clustering,  for  example).  The  variation  of  the  grid 
excludes  the  possibility  of  simply  calculating  the  Taylor  series  coefficients  or  the 
GLA  coefficients  in  the  baseline  design’s  grid,  and  applying  them  to  all  steps  in  an 
iteration  in  the  optimization  process.  This  problem  had  been  previously  identified 
(Ref.  3),  but  was  solved  by  “working  around”  the  problem,  by  changing  the  objective 
function.  This  research  has  identified  a  general  solution  to  this  problem,  by  creating 
a  set  of  “correction  point”  grids,  which  are  independent  of  the  size  of  the  CFD  grid, 
A  system  of  correction  points  was  defined  for  each  wall  surface  (the  inlet  com¬ 
pression  ramp,  for  instance),  in  non-dimensionalized  coordinates,  x/L,  where  L  is  the 
length  of  that  surface.  (See  Fig.  2.2.)  100  to  3500  correction  points  were  used  for 
each  surface,  depending  on  the  length  of  the  surface.  At  the  beginning  of  each  itera¬ 
tion,  the  approximation  concept  coefficients  (e.g.  P(X'’)  and  V  P(X°) )  were  interpo¬ 
lated  from  the  baseline  results,  in  the  CFD  grid  system,  into  the  non-dimensionalized 
coordinates.  As  the  search  progressed  and  the  grid  changed,  £q.  (2.3)  was  applied  at 
the  correction  points  to  find  the  correction  factors.  Wall  surface  pressures  were  then 
calculated  at  these  points,  which  were  then  scaled  to  the  proper  physical  positions. 
No  conversion  back  into  the  CFD  grid  system  was  done.  In  fact,  to  do  so  introduces 
interpolation  and  roundoff  errors  which  were  found  to  yield  incorrect  search  direc¬ 
tions. 

Another  important  consideration  when  using  approximation  concepts  in  aerody¬ 
namic  optimizations  is  that  the  nature  of  the  flow  in  a  given  region  may  strongly  affect 
the  behavior  of  the  correction  factors.  For  example,  the  presence  of  shock  impinge¬ 
ment  points  which  move  as  the  design  is  changed  would  introduce  severe  nonlineari¬ 
ties  into  the  problem  which  such  relatively  simple  approximation  concepts  could  not 
handle.  This  problem  was  encountered  in  optimizing  the  forebody/inlet  region.  (See 
Chapter  5.) 


(a)  Constructing  the  correction  point  grid  system  for  the 
forebody. 


(b)  Converting  the  forebody  wall  pressure  distribution  into 
correction  point  grid  system. 


Fig.  2.2  Example  of  correction  point  grid  system. 
(Only  forebody  shown.) 
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2.5  Summaiy 

Local  approximation  concepts  (based  on  Taylor  series  expansions  of  point  prop¬ 
erties)  and  GLA’s  (which  combine  results  from  detailed  and  approximate  models) 
were  investigated  in  this  research.  Both  methods  seek  to  reduce  the  number  of 
detailed  analyses  required  in  an  optimization.  But  there  are  special  problems  inherent 
in  applying  these  approximation  concepts  to  aerodynamic  problems.  This  research 
has  yielded  a  general  solution  to  the  problem  of  applying  correction  factors  to  grids 
which  change  size  and/or  clustering  as  the  design  is  changed. 


Chapter  3 


Details  of  Flow  Solvers 


3.1  Introduction 

The  analysis  programs  provided  objective  function  and  constraint  data  to  the 
optimizer.  Two  different  levels  of  sophistication  were  used  for  these  programs:  a 
“detailed”  level  and  an  “approximate”  level.  These  levels  of  modeling  can  be  used 
individually  in  the  optimizer,  or  in  concert,  to  form  “Global-Local  Approximations”, 
as  detailed  in  Section  2.3.  The  detailed  analysis  program,  which  is  herein  considered 
the  most  accurate,  consisted  of  a  CFD  flow  solver  for  the  scram  jet  vehicle,  incorporat¬ 
ing  a  1-D,  constant  Mach  number  combustor  model.  Separate  approximate  analysis 
programs  were  used  for  the  nose/forebody  region  and  for  the  nozzle/afterbody  region. 
These  used  the  same  combustor  model  as  the  CFD  program.  The  forebody/inlet 
approximate  flow  solver  was  based  on  oblique  shock  theory,  and  allowed  for  arbitrary 
placement  of  the  shocks  there.  A  simple  Method  of  Characteristics  (MOC)  model  of 
the  nozzle/afterbody  was  used  as  the  approximate  model  of  that  region.  In  addition,  a 
relatively  crude  nozzle/afterbody  model,  using  1 -dimensional  isentropic  flow  theory, 
was  used  as  an  alternative  model,  to  examine  the  effects  of  various  levels  of  fidelity 
in  the  approximate  modeling. 

This  chapter  outlines  the  details  of  each  of  these  flow  solvers,  and  examines  its 
fidelity. 

3.2  CFD  Solver 

The  detailed  analysis  program  consisted  of  a  2-D  inviscid  CFD  solver  for  the 
nose,  forebody  and  inlet,  a  l-D  constant  Mach  number  combustor  model,  and  a  2-D 
inviscid  CFD  solver  for  the  nozzle/aftcrbody  region.  All  flow  in  the  forebody/inlet, 
and  nozzle/afterbody  regions  was  treated  as  being  ideal-gas  frozen  flow,  with  y  -  1.4. 

The  solver  used  the  Steger- Warming  implicit  time  marching  algorithm  for  solv¬ 
ing  the  unsteady  Euler  equations^to  a  steady  state.  (Ref.  26)  This  algorithm  is  an 
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upwind-differenced  fiux-splitting  technique,  using  exact  flux  Jacobians.  As  imple¬ 
mented  in  this  research  effort,  the  scheme  is  second-order  accurate  spatially  every¬ 
where  except  in  the  vicinity  of  large  pressure  gradients  (such  as  near  shocks),  where  a 
“switch”  lowers  the  spatial  differencing  to  first-order.  Generalized  curvilinear  coor¬ 
dinates  were  used  to  allow  various  types  of  geometries  to  be  analyzed  with  the  same 
flow  solver. 

This  well-proven  scheme  has  one  major  drawback,  in  that  it  is  inherently  more 
dissipative  than  some  other  schemes,  such  as  Harten’s  Scheme.  One  result  of  this  is 
that  it  does  not  capture  shocks  as  crisply  as  those  other  schemes.  Thii.  fact  can  have 
an  impact  on  the  optimization  process,  if  a  Global-Local  Approximation  (GLA)  is 
used.  In  such  a  case,  the  detailed  and  approximate  models  may  not  represent  shocks 
equally  crisply,  which  can  lead  to  problems  in  calculating  the  correction  factors  in  the 
GLA  as  shocks  move.  (See  Section  5.5.1.)  Shock  diffusion  caused  no  notijeable 
problems  with  optimizations  not  involving  GLA’s. 

The  flow  solution  process  involved  solving  for  the  flow  in  the  forebody/inlet 
region,  without  regard  to  the  combustor  or  nozzle  flows,  since  supersonic  flow  was 
assumed  throughout.  Then  the  combustor  flow  was  calculated,  using  the  1-D  constant 
Mach  number  model  described  in  Section  3.3.  Finally,  the  nozzle  region  solution  was 
found,  using  the  combustor  exit  flow  and  the  free-stream  flow  from  the  forebody  solu¬ 
tion  as  upstream  boundary  conditions. 

A  special  grid  generation  package  was  written  for  the  scramjet-powered  vehicle 
problem.  A  sample  grid  made  by  this  package  is  shown  in  Fig.  3.1.  A  simplified 
zonal-type  grid  was  used,  with  one  zone  being  defined  by  the  grid  lines  which  enter 
the  engine  and  the  other  by  those  which  do  not.  Details  of  the  grid  are  shown  in  Fig. 
3.2.  This  figure  shows  the  generalized  coordinate  system.  (^,  q).  Note  that  there  is 
elliptically-generated  q -clustering  about  the  walls.  This  peculiarity  for  an  inviscid 
problem  is  brought  about  by  the  need  to  insure  that  the  grid  lines  near  the  wall  are 
very  nearly  tangent  to  the  wall,  when  calculating  the  wall  boundary  conditions.  For 
the  large  flow  contraction  from  the  nose  to  the  inlet,  an  elliptic  grid-generation 
scheme  for  the  t)  coordinates  was  found  to  be  an  effective  way  of  insuring  the  grid 
lines  near  the  wall  were  parallel  to  the  wall,  while  keeping  the  changes  of  slope 
smooth.  This  smoothness  was  found  to  prevent  pressure  “spikes”  near  the  wall  cor¬ 
ners,  such  as  at  the  nose  and  compression  ramp  comers.  Such  spikes  caused  prob¬ 
lems  with  the  optimization  process,  since  they  often  introduced  spurious  changes  in 
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the  objective  function  and  constraints.  Indeed,  for  high  Mach  numbers,  sharp  changes 
in  grid  line  slope  made  the  CFD  solution  unstable. 

Note  the  clustering  of  horizontal  grid  lines  in  the  freestream  ahead  of  the  inlet  in 
Fig.  3.1.  This  is  an  artifact  of  the  simple  zonal  nature  of  the  grid:  the  grid  lines  which 
have  not  yet  entered  the  inlet  are  clustered  in  the  same  way  as  those  which  have 
entered  the  inlet. 


Fig.  3.1  Typical  CFD  grid  for  the  2-D  scramjet 
vehicle. 


The  lines  of  constant  ^  (i.e.  vertical  lines)  were  generated  algebraically,  using 
cosine-type  clustering  near  the  nose,  inlet  compression  ramp,  inlet,  nozzle  throat  and 
nozzle  cowl  lip. 

A  maximum  grid  size  of  400x50  was  found  to  give  a  reasonable  compromise 
between  accuracy  of  the  CFO  solution  and  computational  time.  Actual  grid  sizes 
were  somewhat  smaller  in  the  ^  direction,  depending  on  the  values  of  the  design  vari¬ 
ables.  'Fhe  grid  size  in  the  i\  direction  was  Icept  constant  at  SO  grid  points.  Run  times 
on  an  IBM  RS-6000/320  Powerstation  were  about  2.5  hours. 
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(a)  Forebody/inlct  region. 


10 


Nozzle 


0 


10 


so 


(b)  Nozzle/afterbody  region. 

Fig.  3.2  E>etiiils  of  a  typical  CFD  grid  for  the  2-D 
scrarojet  vehicle. 
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As  the  CFD  solution  developed,  the  nose  and  ramp  shocks  formed  near  the  walls, 
and  then  rotated  to  their  correct  angles  after  further  iteration.  This  progression  is 
shown  in  Fig.  3.3.  During  the  period  when  the  shocks  impinged  upon  the  cowl,  and 
continued  their  progress  (if  the  design  warranted  further  shock  rotation),  the  algorithm 
was  often  found  to  be  unstable.  This  problem  was  circumvented  by  using  the  switch 
which  controlled  the  spatial  accuracy  of  the  scheme  to  impose  first-order  accuracy 
throughout  the  domain  until  the  shocks  had  n:oved  approximately  to  their  final  posi¬ 
tions.  Then  the  switch  was  changed  so  that  the  scheme  was  second-order  accurate 
everywhere  except  very  near  the  shocks. 

A  representative  CFD  solution  is  shown  in  Fig.  3.4.  Note  the  diffusion  of  the 
shocks  at  the  nose,  inlet  compression  ramp  comer,  and  off  the  inlet  cowl  lip.  Recall 
that  the  combustor  is  modeled  as  a  1-D  constant  Mach  number  process,  as  described 
in  Section  3.3.  Convergence  histories  of  the  forebody/inlet  and  nozzle/afterbody  solu¬ 
tions  are  shown  in  Fig.  3.5  and  Fig.  3.6,  respectively. 

Two  measures  of  convergence  were  used:  the  sum  of  the  L-2  norms  of  the 
changes  in  flew  variables  from  one  iteration  to  the  next;  and  the  variation  of  thrust  or 
drag,  as  appropriate  for  the  forebody  or  nozzle,  from  one  iteration  to  the  next.  There 
is  substantial  pre-shock  oscillation  in  the  solutions,  which  keeps  the  L-2  norms  from 
decreasing  more  than  about  1.5  orders-of-  magnitude  for  the  forebody/inlet  solution. 
(See  Fig.  3.5a.)  However,  the  effect  of  this  on  the  drag  results  is  small,  as  shown  in 
Fig.  3.5b.  For  the  design  shown,  after  700  iterations,  the  peak-to-peak  variation  in 
drag  is  less  than  0.05%  of  the  mean  value.  This  level  of  “noise”  can  be  important, 
though,  when  using  the  CFD  solver  to  calculate  gradients  by  a  “brute  force”  finite  dif¬ 
ference  method.  In  such  calculations,  we  perturb  the  design  variables  by  a  small 
amount,  and  difference  the  CFD  results.  To  keep  the  resulting  differences  above  the 
noise  level,  it  was  found  that  the  design  variables  had  to  be  perturbed  by  at  least  1%. 

On  the  other  hand,  the  nozzle  solution  converged  to  well  over  6  orders-of-magni- 
tude  decrease  in  the  L-2  norms;  after  800  iterations,  the  variation  of  thrust  was  negli¬ 
gible.  See  Fig.  3.6a  and  b. 

Based  upon  such  observations,  all  CFD  runs  were  limited  to  a  maximum  of  850 
iterations  for  the  nose/forebody  region,  and  1100  iterations  (or  5  orders-of-magnitude 
decrease  in  L-2  norms)  for  the  nozzle/afierbody  region.  It  is  possible  to  achieve  better 
convergence  with  this  algorithm,  at  the  expense  of  lowering  the  scheme  to  first-order 
accurate  over  more  of  the  flow  domain.  However,  it  was  found  that  doing  this  intro- 
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(b)  Inlet  region 


Hg.  3  A  Mach  number  contours  from  a  representative  CFD 
solution;  (a)  entire  vehicle,  (b)  inlet  region,  and  (c)  nozzle 
region.  ' 
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Fig.  3.4  (cont’d)  CFD  solution  for  optimum 
configuration,  nozzle  region. 


duces  errors  in  the  shock  angles,  which  can  have  a  dramatic  effect  on  the  thrust  and 
drag  calculations,  due  to  the  sensitivity  of  the  results  to  shock  reflections  near  the 
inlet. 

To  validate  the  CFD  solver,  the  nozzle/afterbody  configuration  of  Ref.  4  (p.  723, 
Fig.  4.)  was  modeled,  and  the  results  were  compared  to  the  experimental  results 
quoted  in  that  reference.  A  comparison  of  the  pressure  distributions  along  the  noz¬ 
zle  ramp,  for  the  Steger-Warming  CFD  code,  and  for  the  experimental  results,  is 
shown  in  Fig.  3.7.  The  results  are  generally  quite  good,  especially  bearing  in  mind 
that  the  CFD  code  is  for  inviscid  flow. 

Validation  of  the  CFD  solver  for  the  nose  and  forebody  was  done  by  comparing 
the  results  to  those  from  oblique  shock  theory.  An  example  of  such  a  comparison  is 
shown  in  Fig.  3.8.  Again,  the  agreement  is  quite  good.  But  note  the  diffusion  of  the 
^ocks,  due  both  to  the  discretization  of  the  CFD  method  and  to  the  dissipative  nature 
of  the  Steger-Warming  algorithnu 
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Fig.  3.5  CFD  convergence  histoiies  for  foiebody/inlet  region 
for  a  typical  design;  (a)  sum  of  L-2  norms,  and  (b)  drag. 
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(a)  History  of  nozzle/afterbody  L-2  norms. 


Fig.  3.6  CFD  convergence  histories  for  nozzlc/afterbody  region 
for  a  typical  design;  (a)  Sum  of  L>2  norms,  and  (b)  thrust. 
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Fig.  3.8  Comparison  of  CFD  nose/forebody  wall 
pressure  distribution  with  oblique  shock  theory. 


3.3  Combustor  Model 

The  combustor  was  modeled  as  a  1 -dimensional,  constant-Mach  number  process. 
The  fuel-air  mixture  was  assumed  to  be  a  perfect  gas  with  y  =  i  •4.  The  addition  of 
mass  due  to  fuel  injection  was  neglected.  It  would  be  more  customary  to  employ  a 
constant-pressure  combustor,  but  as  Ref.  22  points  out,  the  constant-Mach  number 
combustor  offers  two  advantages  from  an  analysis  standpoint:  it  explicitly  avoids 
choking  in  the  combustor,  and  it  greatly  simplifies  the  analysis.  In  fact,  the  differen¬ 
tial  equations  for  generalized  steady  1 -dimensional  flow  with  heat  addition,  as  given 
in  Ref.  32,  can  be  reduced  to  algebraic  equations  for  constant-Mach  number  flow. 
Note  that,  to  maintain  a  constant  Mach  number  while  heat  is  being  added  to  the  flow, 
the  combustor  cross-sectional  area  must  expand  from  inlet  to  exit.  Such  a  combustor 
would  be  difficult  to  employ  in  practice,  as  it  would  require  that  the  heat  addition  be 
specifically  tailored  along  the  length  of  the  combustor.  However,  for  the  purposes  of 
demonstrating  different  optimization  strategies,  the  constant-Mach  number  combustor 
provides  a  useful,  simple  tool  which  captures  the  salient  features  of  total  temperature 
increase  and  total  pressure  decrease. 

The  combustor  process  is  governed  by  the  following  relations,  in  which  the  sub¬ 
scripts  “i”  and  “e”  refer  to  conditions  at  the  combustor  inlet  and  exit,  respectively. 


Y  A  2/(l*yM]) 
—  =  (— ) 


Pi 


Ai  T, 
—  (—) 
A,T/ 


1/2 


(3.1) 

(3.2) 


where  T,  and  p,  are  the  total  temperature  and  total  pressure,  A  is  the  combustor 
cross-sectional  area,  and  is  the  (constant)  combustor  Mach  number. 
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The  area  ratio  was  set  by  determining  a  baseline  value  of  total  temperature  ratio 
forgiven  freestream  and  combustor  conditions,  and  applying  (3.1)  along  with  the  fol¬ 
lowing  relations  for  r, . 

T.  »  r,  +  7^  (3.3) 


where 


Q  =  (loicA^mu  (3.4) 

in  which  Q  is  the  heat  added  by  combustion,  is  the  heating  value  of  the  fuel,  (|>  is 
the  equivalence  ratio,  is  the  stoichiometric  fuel/air  mass  ratio,  and  is 

the  mixing  efficiency.  The  freestream  conditions  were  set  equal  to  those  for  A/.  =  6.0 
at  30  km  altitude,  and  the  combustor  Mach  number  was  set  at  1.20.  The  fuel  was  gas¬ 
eous  hydrogen,  and  was  set  at  0.90.  ITie  resulting  area  ratio  was  about  2.693;  for 
simplicity,  and  because  the  choice  is  actually  somewhat  arbitrary,  it  was  set  to  2.5  for 
all  optimizations. 

Note  that  the  area  expansion  in  the  combustor  results  in  a  large  portion  of  the 
thrust  coming  from  the  combustor;  for  most  designs  evaluated  in  the  optimizations, 
“combustor  thniist”  accounted  for  about  50%  of  the  total  thrust. 

Although  the  combustor  length  is  not  critical  to  the  design  optimization  for  the 
problem  at  hand,  it  was  felt  that  a  reasonable  length  should  be  used.  The  combustor 
length  was  chosen  by  specifying  that,  for  all  designs  used  in  the  optimizations,  the 
combustor  should  be  at  least  long  enough  to  guarantee  that  mixing  and  complete  com¬ 
bustion  take  place  within  the  combustor.  Using  H2-air  reaction  rate  data  from  Ref.  35, 
mixing  efficiency  data  from  Ref.  19,  and  results  from  the  initial  design  used  in  the 
optimizations,  the  minimum  combustor  length  was  found  to  be  about  4  meters.  To 
guarantee  that  all  designs  that  would  be  evaluated  in  the  optimizations  would  have 
complete  combustion,  a  combustor  length  of  10  meters  was  used.  The  combustor 
performance  is  characterized  in  Fig.  3.9.  Shown  in  Fig.  3.9a  are  curves  of  total  tem¬ 
perature  ratio  and  total  pressure  ratio  across  the  combustor,  as  a  function  of  combustor 
Mach  number.  Fig.  3.9h  is  a  parameterization  of  combustor  Mach  number  for  a  range 
of  ramp  angles,  assuming  a  constant  8.0  degree  nose  angle.  Perfectly-placed  fore¬ 
body/inlet  shocks  were  assumed  here.  (This  restriction  was  removed  for  the  actual 
optimizations.)  Fig.  3.9c  is  a  similar  parameterization,  for  the  thrust  produced  by  the 


Fig.  3.9  Combustor  performance:  (a)  combustor  total 
temperature  and  total  pressure  ratios  vs  combustor 
Mach  number;  (b)  combustor  Mach  number  vs  ramp 
angle,  for  constant  nose  angle:  and  (c)  combustor  force 
vs  ramp  angle  foi^  constant  nose  angle. 


combustor  walls.  It  is  obvious  that  there  is  a  great  deal  to  be  gained  (as  far  as  com¬ 
bustor  performance  is  concerned)  by  using  a  large  degree  of  flow  turning  in  the  fore¬ 
body.  For  the  nose/forebody  optimizations  the  nose  angle  was  limited  to  a  maximum 
of  11.0  deg,  and  the  ramp  angle  was  limited  to  20.0  deg,  since  improper  placement  of 
the  cowl  shock  (which  is  allowed  by  the  optimization  strategy  chosen  for  this 
research)  could  result  in  engine  unstart  for  angles  much  greater  than  these.  As  can  be 
seen  from  Fig.  3.9,  this  greatly  limits  the  combustor  performance  that  can  be  achieved 
by  the  optimizations.  However,  allowing  for  arbitrary  shock  placement  helps  us 
examine  more  thoroughly  the  performance  of  the  optimization  techniques  in  various 
flow  regimes.  The  limitation  in  combustor  performance  is  thus  justified. 


3.4  Approximate  Forebody/Inlet  Model 

The  ’“approximate”  method  used  to  model  the  nose,  compression  ramp  and  inlet 
was  based  on  oblique  shock  theory.  The  conventional  design  process  for  such  a  con¬ 
figuration  would  limit  the  possible  designs  to  those  with  perfect  shock  placement, 
with  the  nose  and  compression  ramp  shocks  hitting  the  cowl  lip,  and  the  cowl  shock 
impinging  on  the  inlet  comer.  But  in  this  demonstration  of  advanced  optimization 
techniques,  the  optimization  process  proceeds  without  explicit  regard  to  the  shock 
locations.  Therefore,  dii.s  approximate  model  had  to  be  given  the  capability  to  treat 
designs  which  did  not  perfectly  place  the  shocks. 

This  approximate  model  begins  by  calculating  the  shock  strengths  and  angles  for 
the  nose  and  compression  ramp  shocks  (see  Fig.  3.10.),  using  oblique  shock  theory. 
It  then  determines  into  which  of  four  general  categories  the  shock  structure  will  fall: 

Case  1.  Nose  shock  and  ramp  shock  do  not  intersect  each 

other  before  reaching  the  cowl  lip,  and  neither  shock 
impinges  anywhere  on  the  cowl.  See  Fig.  3.10a. 

Case  2.  Nose  shock  and  ramp  shock  do  not  intersect  before 
reaching  the  cowl  lip,  and  both  hit  the  cowl  surface 
either  af  or  downstream  of  the  lip.  See  Fig.  3.10b. 
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Case  3.  Nose  shock  and  ramp  shock  do  not  intersect  before 
reaching  the  cowl  lip,  and  the  nose  shock  does  not 
intersect  the  cowl,  while  the  ramp  shock  hits  the  cowl 
either  at  or  downstream  of  the  Up.  See  Fig.  3.10c. 

Case  4.  Nose  shock  and  ramp  sliock  intersect  before  reaching 
the  cowl  lip.  See  Fig.  3.10d. 

For  Case  1,  the  model  then  calculates  the  strength  and  angle  of  the  cowl  shock.  It 
finds  the  intersection  of  this  shock  with  the  fuselage  surface,  and  determines  if  the 
shock  will  reflect  from  the  fuselage  surface  before  passing  into  the  inlet  plane.  If  so, 
the  reflection  properties  are  calculated,  and  the  intersection  of  the  reflected  shock 
with  the  inlet  plane  is  found.  If  not,  the  intersection  of  the  cowl  shock  with  the  inlet 
plane  is  found.  The  model  then  calculates  surface  pressures  along  the  fuselage  fore< 
body  and  cowl  Up,  and  finds  the  flow  properties  along  the  inlet  plane.  This  informa¬ 
tion  is  used  to  calculate  forebody  and  cowl  drag  and  the  average  inlet  plane  flow 
properties,  which  are  passed  to  the  combustor  model. 

For  Case  2,  no  shock  emanates  from  the  top  of  the  cowl  Up,  since  the  flow  there 
is  already  parallel  to  that  surface.  Thus  the  intersection  of  the  nose  shock  with  the 
cowl  surface  is  found,  and  the  reflection  is  found.  This  reflection  then  intersects  the 
ramp  shock.  The  method  outlined  in  Ref.  29  for  finding  the  intersection  of  two 
shocks  of  opposite  families  is  used  to  find  the  resulting  two  shocks  and  slip  line.  The 
model  then  determines  if  these  shocks  and  slip  line  intersect  the  inlet  plane  directly,  or 
reflect  off  the  inlet  surfaces  first.  It  calculates  such  reflections,  if  necessary,  and  so 
finds  the  average  inlet  plane  flow  pmoerties  and  surface  pressure  distributions. 

In  Case  3,  there  is  a  shock  off  the  cowl  Up.  The  model  finds  the  intersection  of 
this  shock  with  the  ramp  shock,  and  finds  any  subsequent  intersections  of  shocks  with 
the  inlet  surfaces  and  then  wi^h  the  inlet  plane.  It  then  determines  the  surface  pres¬ 
sures  and  inlet  plane  flow  properties. 

Note  that,  in  Cases  1-3,  the  model  treats  the  intersection  of  slip  lines  with  shocks 
very  crudely.  It  simply  terminates  a  slip  line  when  it  reaches  a  shock.  Also,  the  effect 
of  the  slip  line  on  the  strength,  shape  and  direction  of  downstream  shocks  is 
neglected.  The  only  real  accounting  of  a  slip  line  occurs  when  intersects  the  inlet 
plane.  In  such  a  case,  the  difference  in  flow  properties  on  either  side  of  the  slip  line 
affect  the  average  inlet  plane  flow  properties  that  are  passed  to  the  combustor  model. 
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Fig.  3.10  Forebody/inlet  shock  configurations  used  u) 
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Case  4  invokes  the  calculation  of  the  intersection  of  two  shocks  of  the  same  fam¬ 
ily,  as  explained  in  Ref.  29.  It  determines  the  resulting  “transmitted”  and  “reflected” 
waves  and  slip  line.  Then  it  checks  to  see  if  the  transmitted  shock  hits  the  cowl. 
From  this  point,  it  proceeds  in  much  the  same  manner  as  Cases  1-3.  Note  that  the  slip 
line  from  the  intersection  of  the  nose  and  compression  ramp  shocks  is  tracked  all  the 
way  to  the  inlet  plane.  Furthermore,  the  treatment  of  the  intersection  of  this  slip  line 
with  a  shock  is  rather  better  than  similar  intersections  in  Cases  1-3,  in  that,  upon 
crossing  the  slip  line,  the  shock  is  deflected.  This  additional  fidelity  (and  complexity) 
was  warranted  by  the  influence  tliis  slip  line  has  on  the  location  of  the  intersection  of 
the  cowl  shock  with  the  wall  surface.  This  intersection  has  a  great  effect  on  the  drag 
and  thrust  results,  since  the  presence  of  a  shock  reflection  near  the  inlet  corner  can 
greatly  increase  fombody  drag,  as  well  as  the  average  pressure  and  density  in  the  inlet 
plane. 

Note  also  that  the  “reflected”  or  “secondary”  wave  from  the  shock  intersection  in 
Case  4  is  tracked  no  further  than  its  intersection  with  the  wall.  This  was  done  to  sim¬ 
plify  the  analysis. 

Since  this  approximate  model  uses  oblique  shock  theory  for  its  shock  calcula¬ 
tions,  the  shock-j'inip  values  and  shock  angles  for  the  nose,  cormpression  ramp  and 
cowl  shocks  may  actually  be  slightly  more  accurate  than  their  CFD  counterparts. 
However,  the  approximations  concerning  the  treatment  of  slip  lines  in  Cases  1-3  and 
the  treatment  of  the  reflected  wave  in  Case  4  introduce  errors  in  the  results  of  the 
approximate  model. 

The  greatest  source  of  error  in  the  approximate  forebody  model  is  that  it  does  not 
model  the  interaction  of  the  inlet  shock  system  with  the  expansion  at  the  inlet  comer. 
It  simply  assumes  that  the  effect  of  the  expansion  is  to  “freeze”  and  to  average  the 
inlet  plane  flow  conditions.  This  assumption  greatly  simplifies  the  model.  Compari¬ 
sons  of  the  approximate  model’s  results  with  those  of  the  CFD  solver  (which  models 
the  flow  through  the  entire  inlet  length)  show  that  this  assumption  is  reasonable,  if 
somewhat  in  error.  Ip  the  optimizations,  the  largest  difference  between  drag  values 
observed  for  the  two  different  models  was  less  than  12%,  and  the  largest  difference 
between  thrust  values  was  less  than  23%.  Most  often,  the  differences  were  less  than 
7%.  The  values  of  net  thrust,  however,  were  in  error  by  as  much  as  80%.  But  for  the 
purposes  of  optimization,  this  model  docs  quite  well,  as  will  be  shown  in  Section  5.5. 
It  correctly  captures  the  trends  inmet  thrust  that  changes  in  the  design  variables  cause. 
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This  was  demonstrated  by  the  similar  results  obtained  from  optimizations  using  CFD 
alone  and  using  the  approximate  model  alone. 

3.5  Approximate  Models  for  the  Nozzle 

For  the  purpose  of  comparing  models  of  differing  accuracy,  two  different  approx¬ 
imate  models  of  the  nozzle/afterbody  were  made.  The  first  was  a  straightforward 
application  of  the  steady,  2-dimensional  Method  of  Characteristics  (MOC),  as  detailed 
in  Ref.  33.  The  second  was  a  simple  application  of  1 -dimensional  isentropic  flow, 
assuming  that  no  shock  waves  impinge  upon  the  nozzle  ramp. 


3J.11  MOC  Model 

A  typical  MOC  mesh  for  the  scramjet  nozzle/afterbody  is  shown  in  Fig.  3.12a. 
In  this  model,  all  characteristics  emanate  from  the  expansion  comers.  That  is,  half  of 
tine  characteristics  have  their  origins  at  the  ramp  comer,  and  the  other  half  have  their 
origins  at  the  cowl  comer.  Each  expansion  fan  is  divided  into  50  characteristic  lines. 
(For  clarity  in  Fig.  3.12a  ,  only  25  characteristic  lines  are  shown  emanating  from  each 
comer.) 

Depending  on  how  much  the  flow  at  the  cowl  lip  is  under-  or  over-expanded, 
there  will  be  either  an  oblique  shock  or  an  expansion  wave  there  to  adjust  the 
freestreara  flow  pressure  and  direction,  and  possibly  also  one  to  adjust  the  nozzle 
flow.  See  Fig.  3.1 1.  It  is  assumed  in  the  MOC  analysis  that  no  such  shock  intersects 
the  nozzle  ramp.  This  is  justified  because,  for  all  the  optimiztions  performed,  the  flow 
was  only  slightly  over-expanded  or  under-expanded.  Also  in  the  MOC  analysis,  the 
interactions  between  the  characteristics  and  these  waves  emanating  from  the  cowl  lip 
are  neglected.  This  is  justified  on  the  assumption  that  the  resulting  reflection  of  the 
characteristic  is  unlikely  to  intersect  the  ramp,  given  out  geometry.  (See  Fig.  3.12a.) 
Also,  as  pointed  out  in  Ref.  29,  Aiticle  16.8,  the  result  of  an  intersection  between  an 
expansion  wave  and  a  shock  of  the  same  family  is  a  weak  reflected  expansion  wave 
and  a  curving  of  the  shock.  The  weak  expansion  wave  is  not  likely  to  have  much  of  an 
impact  on  the  ramp  surface  pressure  distribution,  even  if  it  is  felt  there. 

Fig.  3.12a  also  shows  that,  since  the  characteristics  emanate  only  from  the  cor¬ 
ners,  no  characteristics  intersect  the  ramp  between  the  ramp  comer  and  the  intersec- 
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tion  of  the  first  reflected  characteristic  with  the  ramp.  The  correct  wall  pressure  will 
be  found  at  the  first  intersection  with  the  ramp,  but  the  variation  of  pressure  between 
the  ramp  corner  and  this  point  will  be  unknown.  In  this  model,  the  wall  pressure  is 
assumed  to  vary  linearly  between  these  two  pointe. 

Fig.  3.12b  compares  wall  pressure  distributions  for  the  MOC  model  with  those 
from  the  CFD  solver,  for  a  representative  configuration.  Note  the  linear  distribution 
of  pressure  immediately  downsueam  of  the  ramp  comer,  as  mentioned  in  the  previous 
paragraph.  This  is  due  to  the  lack  of  characteristics  modeled  in  this  region.  Note,  also 
that  the  general  character  of  the  CFl)  wall  pressure  distribution  is  reasonably  well  rep¬ 
licated  by  the  MOC. 

Although  the  values  of  net  thmst  found  by  the  MOC  model  differ  from  the  CFD 
values,  the  general  trends  of  behavior,  as  the  design  variables  are  allowed  to  vary,  do 
follow  those  of  the  CFD  quite  well,  as  shown  in  Fig.  4.2  and  Fig.  4.3. 
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Fig.  3.11  Example  of  waves  emanating  from  cowl  lip 


(a)  MOC  Nozzie  Mesh. 


(b)  Ramp  Pressure  Distribution. 

Fig.  3.12  Method  of  Characteristics  nozzle  mesh 
and  ramp  wall  pressure  distribution,  with 
comparison  to  CFD  results.  (Same  configuration  as 
in  Rg.  3.4.) 
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3.5.2  l-D  Isentropic  Flow  Model 

The  1-D  flow  model  assumes  the  flow  in  the  nozzle/aftcrbody  region  behaves 
one-dimensionally,  with  flow  cross-sectional  areas  defined  by  the  area  between  the 
ramp  surface  and  either  the  cowl  inner  surface  (if  appropriate),  or  the  “contact  sur¬ 
face”  bounding  the  nozzle  plume.  See  Fig.  3.13a. 

This  contact  surface  is  found  by  examining  the  CFD  solution  for  a  baseline 
design  configuration.  To  identify  the  location  of  the  contact  surface  at  some  station 
downstream  of  the  cowl  lip,  we  draw  a  vertical  line  segment  at  that  station  and  search 
along  this  line,  between  the  lowest  point  in  the  freestream  and  the  point  which  is  on  a 
level  with  the  cowl  Up.  See  Fig.  3.11.  We  identify  the  point  of  maximum  fluid  den¬ 
sity  as  the  point  on  the  contact  surface.  This  definition  of  contact  surface  is  motivated 
by  the  observation  that,  for  the  configurations  used  in  the  optimizations,  the  two 
waves  emanating  from  the  cowl  lip  are  usually  both  weak  shocks.  We  know  that  the 
contact  surface  must  lie  somewhere  between  these  two  waves,  and  that  the  point  of 
maximum  density  will  also  lie  between  them.  Since  this  point  is  readily  identifiable 
from  the  CFD  solution,  the  point  of  maximum  density  is  attractive  for  identifying  the 
contact  surface. 

When  the  contact  surface  is  found  by  this  method,  it  has  a  “staircase”  quality,  as 
shown  in  Fig.  3.13a.  This  is  due  to  the  discrete  nature  of  the  CFD  solution.  For  the 
1-D  model,  this  surface  is  smoothed  by  connecting  the  midpoints  of  the  vertical  legs 
of  the  steps  as  shown. 

A  typical  nozzle  ramp  pressure  distribution  found  by  the  1-D  isentropic  flow 
model  is  shown  in  Fig.  3.13b.  The  CFD  solution  is  also  shown.  Note  that  over  the 
first  half  of  the  length,  the  1-D  model  does  not  provide  a  good  representation  of  the 
flow  at  the  wall,  as  would  be  expected  from  its  inherent  limitations.  However,  over 
the  downstream  half  it  becomes  much  more  accurate. 

Reference  to  Fig.  4,3  and  Fig.  4.4  show  that,  as  expected,  the  net  thrust  predic¬ 
tions  from  the  1-D  model  do  not  follow  the  CFD  values  as  well  as  the  MOC  predic¬ 
tions  do. 
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Oist.  Downstream  of  Throat, 
(a)  Contact  surface. 


(b)  1-D  ramp  pressure  distribution. 


Fig.  3.13  Typical  results  from  1-D  isentropic  flow 
model:  (a)  contact  surface,  and  (b)  ramp  surface 
pressure  distribution,  compared  to  C^. 
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3.6  Summary 

A  2-D  Euler-solver  CFD  program  provided  the  “detailed”  analysis  capability  for 
the  optimizations.  The  “approximate”  model  of  the  nose/forebody  consisted  of  an 
oblique  shock  theory  model  which  allowed  for  simplified  treatment  of  arbitrary  shock 
configurations.  Two  approximate  models  of  the  nozzle/afterbody  region  were  made, 
for  purposes  of  comparison:  one  based  on  the  steady  2-D  Method  of  Characteristics, 
and  a  more  crude  one  based  on  1-D  isentropic  flow  theory.  The  approximate  models 
were  all  shown  to  provide  generally  reasonable  fidelity  in  relation  to  ;:he  CFD  pro¬ 
gram. 

The  detailed  and  approximate  models  were  mtended  to  be  used  individually  in 
the  optimizations,  as  well  as  in  concert,  to  form  the  bases  of  the  “Global-Local 
Approximations”  (GLA’s). 


/ 


Chapter  4 


Noz2ie/Afterbody  Optimizations 


4. 1  Introduction 

The  approximation  concepts  were  first  applied  to  optimizations  of  the  nozzle/ 
afterbody  region.  This  region  is  characterized  by  wall  flow  that  is  essentially  isen- 
tropic,  and  thus  is  a  rather  benign  environment  for  the  approximation  concepts.  Noz¬ 
zle  optimizations  are  also  cheaper  to  perform  than  forebody/inlet  optimizations,  since 
the  nozzle  inlet  conditions  can  be  specified,  and  only  the  nozzle  portion  of  each  of  the 
flow  solvers  has  to  be  run.  Thus,  it  is  particularly  attractive  to  begin  the  optimizations 
with  the  nozzle. 

This  chapter  describes  the  nozzle  optimization  setup.  Then  it  investigates  the 
fidelity  of  each  of  the  approximation  concepts  by  performing  a  series  of  single- vari¬ 
able  design  parameterizations.  Finally,  it  discusses  the  actual  optimizations  that  were 
performed  on  the  nozzle  using  each  of  the  approximation  concepts,  as  well  as  CFD- 
alone  and  each  of  the  uncorrected  approximate  methods.  Forebody/inlet  optimiza¬ 
tions  will  be  treated  in  Chapter  5. 

4.2  General  Problem  Considerations 

Fig.  4. 1  shows  the  setup  for  this  design  problem.  Note  that  the  ramp  surface  is 
defined  by  a  quadratic  equation  in  x,  where  x  is  the  horizontal  distance  from  the 
throat.  (This  follows  from  the  observation  in  Ref.  4  that  skc\^'ed  parabolic  contours 
have  been  considered  for  2-0  NASP  nozzles.  Other  families  of  surface  descriptions 
may  be  just  as  easily  used  in  the  optimization  process.) 

The  following  design  parameters  were  considered  to  have  predominant  influ¬ 
ences  on  the  thrust:  initial  angle  of  the  ramp  (a);  curvature  of  the  ramp,  as 

defirsd  by  the  quadratic  coefficient  in  the  equation  for  the  ramp  surface;  cowl  angle 
O).  Note  that  these  parameters  can  be  expressed  by  combinations  of  other  parame¬ 
ters,  and  that  using  such  alternative  parameters  for  the  design  variables  may  actually 
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result  in  better  optimizer  behavior.  (See  Section  4.4.1)  Other  potentially  important 
design  parameters,  such  as  cowl  length  and  cowl  curvature  were  not  treated,  in  the 
interest  of  limiting  the  size  of  the  problem. 


(fixed) 

Fig.  4.1  Nozzle  optimization  geometry. 


The  approximate  analysis  models  used  were  the  1-D  isentropic  flow  model,  and 
the  Method  of  Characteristics  (MOC)  model.  Local  (Taylor  scries),  and  Global-Local 
Approximations  (GLA’s)  were  used.  (See  Sections  2.2  and  2.3.)  The  CFD  method  is 
herein  deemed  to  be  the  most  accurate  of  the  methods,  since  it  most  accurately  treats 
the  effects  of  the  nozzle  plume,  and  since  it  is  much  mors  finely  disci etized  than  the 
next  most  accurate  method,  the  MOC.  (The  CFD  model’s  accuracy  is  somewhat  less¬ 
ened  by  its  numerical  damping,  which  is  absent  in  the  MOC.)  Thus,  the  CFD  results 
will  be  the  baseline  against  which  all  other  methods  will  be  compared. 
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Throat  flow  conditions  were  fixed,  and  corresponded  to  those  resulting  from  the 
design  shown  in  Fig.  3.4.  (The  forebody/inlet  design  in  that  figure  is  close  to  the  opti¬ 
mum  design,  as  will  be  seen  in  Chapter  5.)  The  throat  flow  was  assumed  to  be  uni¬ 
form,  1-D  flow.  All  flow  in  the  nozzle  was  assumed  to  be  frozen,  perfect  gas  flow, 
with  Y  =  1-4.  The  inflow  conditions  for  the  external  flow  were  set  equal  to  the 
freestream  values. 


4.3  Single-Variable  Parameterizations 

In  order  to  gain  insight  into  the  relative  accuracies  of  the  different  analysis  mod¬ 
el?  and  approximation  concepts,  a  set  of  nozzle  performance  curves  were  generated, 
in  which  the  normalized  net  thrust  is  presented  as  a  function  of  each  design  variable, 
with  the  other  design  variables  held  constant.  Each  of  the  analysis  models  and  approx¬ 
imation  concepts  was  used  and  the  results  were  plotted  on  the  same  plot.  I'hese 
‘‘parameterizations’’  are  shown  in  Fig.  4.2  through  Fig.  4.4. 

In  those  figures,  the  darkened  square  data  points  are  from  the  CFD  solutions. 
The  Global-Local  Approximations  (GLA’s)  using  the  MOC  and  1-D  models  are 
lat>eled  MOC  GLA  and  I  D  GLA.  Both  of  these  GLA’s  and  the  Taylor  series  model 
were  based  upon  baseline  designs  at  a  «  18.0”,  =  18.0”,  and  ramp  curvature  coeffi¬ 

cient  of  0.005.  In  the  plots  the  minor  vertical  displacements  of  the  Taylor  series  and 
GLA  curves  relative  to  the  CFD  baseline  design  solution  are  due  to  the  different  grids 
and  integration  schemes  used  by  these  methods.  These  “DC  shifts’’  have  no  effect  on 
the  trends  of  the  results,  and  thus  do  not  affect  the  optimization  process.  The  shifts  in 
F,,,  at  the  design  point  are  less  than  0.15. 

Figs.  4.2-4.4  highlight  the  advantage  of  the  GLA  approach  over  the  uncorrected 
approximate  methods  and,  to  a  lesser  degree,  over  the  Taylor  series  method.  Near  the 
baseline  design  point,  the  trends  of  the  GLA  results  agree  quite  well  with  those  from 
the  CFD,  even  when  the  uncoirected  approximate  methods  do  not.  Indeed,  in  Fig. 
4.2,  for  a  between  17.5  and  18.5  deg,  the  Taylor  series,  1-D  GLA  and  MOC  GLA 
results  follow  the  CFD  Jesuits  so  well  that  they  are  indistinguishable  from  the  CFD 
curve.  As  the  design  is  shifted  further  away  from  the  baseline,  the  GLA  and  Taylor 
results  still  follow  the  general  trends  of  the  CFD  far  better  than  the  uncoirected 
^proximate  results  do.  Note  in.  Fig.  4,3  that  the  1-D  GLA  and  MOC  GLA  overlay 


Fig.  4.2  Parameterization  of  F  ,  vs  nozzle 
ramp  angle,  a,  for  various  models. 


Fig.  4.3  Parameterization  of  vs  nozzle  ramp 
curvature  coefficient,  ,  for  various  models. 
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Fig.  4.4  Parameterization  of  vs  nozzle  cowl  angle  for 
various  analysis  models. 


each  other  perfectly  over  the  entire  mnge  of  shown. 

The  parameterization  curves  show  that,  for  the  parameterizations  in  this  particu¬ 
lar  problem,  the  GLA's  are  better  at  mimicking  the  CFO  solver  than  the  unconected 
approximate  models  arc.  There  is  really  very  little  to  choose  between  these  two 
GLA’s.  The  Taylor  series  also  does  better  than  the  approximate  methods,  except  in 
the  parameterization,  where  the  uncorrected  MOC  conforms  better  to  the  shape 
of  the  CFD  curve,  for  valiucs  of  far  from  the  baseline  design. 

Note,  especially  in  Fig.  4.3,  that  the  first-order  Taylor  series  approximation 
applied  to  point  properties  (wall  pressures)  does  not  produce  a  linear  behavior. 
Had  we  applied  the  Taylor  series  directly  to  the  integral  property  (F^„),  the  curve 
would  have  been  linear.  This  is  a  clear  example  of  how  nonlinearities  in  a  problem 
can  be  better  huidled  by  some  formulations  of  approximation  concepts  than  by  others. 
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4.4  2~Variable  Nozzle  Optimizations 

Having  established  the  general  behavior  of  the  various  modeling  levels  and 
approximation  concepts,  the  nozzle  design  was  optimized,  using  each  of  them. 
Because  the  effect  of  cowl  angle  on  thrust  was  seen  to  be  much  smaller  than  that  of 
ramp  angle  and  curvature,  only  the  nozzle  ramp  geometry  was  considered  for  deter¬ 
mining  the  design  vaniables. 

Each  optimization  was  performed  by  combining  the  well-proven  optiraiz  ition 
software  package  “Design  Optimization  Tools”,  or  DOT  (See  Ref.  12.),  with  one  or 
more  2-D  flow  solvers  (or  “analysis  programs”).  DOT  u.ses  the  “Modified  Method  of 
Feasible  Directions’"  to  solve  constrained  nonlinear  optimization  problems.  In  this 
case,  DOT  was  instructed  to  find  the  vehicle  design  which  maximized  net  thrust 
(Thrust  -  Drag). 

4.4.1  Cholceof  Design  Variables 

It  has  long  been  known  that  the  choice  of  design  variables  can  have  an  impact  on 
the  optimization  process.  Although  a  and  seem  to  be  the  most  obvious 

choices  for  ramp  design  variables,  another  set  of  design  variables  was  considered. 
This  set  comprised  a  and  A,,,  the  height  of  the  nozzle  exit  plane,  with  being 
determined  from  these  two. 

To  investigate  the  appropriateness  of  each  of  these  2  choices  of  design  variables, 
a  set  of  optimizations  was  run  with  each,  using  the  MOC  only  and  the  conditions 
listed  in  Section  4.4.2.  The  results  arc  shown  in  TABLE  4.1. 

As  can  be  seen  in  TABLE  4.1,  there  was  virtually  no  difference  in  the  optimum 
designs.  However,  it  was  decided  to  use  a  and  h„^  for  all  subsequent  optimizations, 
for  the  following  rea::ons. 

First,  in  order  to  minimize  axial  momentum  loss  in  the  nozzle,  we  would  like  to 
keep  the  nozzle  exit  wall  angle  as  close  to  horizontal  as  possible,  while  still  meeting 
the  other  requirements  of  the  optimization,  such  as  the  fixed  nozzle  length  and  maxi¬ 
mum  nozzle  exit  height.  This  tendency  can  be  incorporated  into  the  optimization 
process  if,  for  a  reasonable  initial  design,  an  increase  in  a  is  accompanied  by  an 
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TABLE  4.1 _ 

Effects,  of  using  different  design  variable  sets  on  the 


nozzle  optimization,  using  uncorrccted  MCKT. 


a  and 

a  and 

Initisi  a 

18.000 

18.000 

IniliJil 

0.0050 

0.0050 

18,051 

18.084 

0,00053 

0.00036 

increase  in  lliat  is.  if  the  effect  of  an  increased  initial  wall  angle  is  tempered 

by  an  increase  in  the  c urvatuie.  (See  Fig.  4.1.) 

To  make  this  design  philosophy  a  part  of  the  optimization  process,  this  tendency 
should  be  reflected  in  the  gradient  calculations,  to  insure  that  a  proper  search  direction 
is  always  chosen.  This  is  the  case  for  the  design  variable  pair  a  and  h„^.  That  is,  if 
we  increase  a  while  Iceeping  constant  (as  is  done  when  calculating  the  partial 
derivative  of  with  respect  to  a)  then  will  necessarily  increase.  We  could 
consider  this  to  be  a  “smart”  design  variable  s<n.  If,  on  the  other  hand,  a  and 
are  used  as  design  variables,  the  two  can  change  independently,  and  this  design  phi- 
losophy  would  be  lost. 

Secondly,  choosing  a  length  instead  of  an  angle  can  help  to  reduce  the  non-lin¬ 
earity  m  the  problem.  As  a  greatly  simplified  example  of  this  effect,  consider  a 
straight  ramp  with  a  constant  pressure  distribution  as  shown  in  Fig.  4.5.  The  axial 
force  exerted  on  the  rarnp  is 


=  pLuma  (4.2) 

which  is  clearly  nonlinear  in  the  design  variable,  a.  However,  if  we  use  as  the 
design  variable,  we  have 


F^r^ph 


(4.3) 
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v/hich  is  linear  in  the  design  variable.  Although  the  actual  nozzle  design  problem  is 
much  more  nonlinear  than  this  example,  we  might  expect  some  improvement  in  the 
optimization  process  due  to  choosing  the  right  design  variables. 


L  - 


F,  =  pttanG  =  ph 


Fig.  4.5  Example  of  the  linearizing  effect  of  choosing 
as  a  design  variable. 


4.4.2  Optimiziition  Setup 

For  each  of  the  nozzle-only  optimizations,  the  design  was  optimized  using  the 
initial  ramp  angle,  a,  and  the  nozzle  exit  height,  as  design  variables.  Optimiza¬ 
tions  were  run  for  all  analysis  models  and  approximation  concepts:  CFD  alone;  1-D 
isentropic  flow  alone;  MOC  alone;  Taylor  series;  GLA  using  1-D  isentropic  flow; 
GLA  using  MOC.  The  resulting  designs  were  compared,  and  the  savings  in  cotnputf  r 
resources  for  each  method  were  then  found.  The  figure  of  merit  for  these  savings  was 
the  number  of  required  calls  to  the  CFD  solver,  exclusive  of  the  calls  used  to  calculate 
the  gradients  of  the  objective  function  and  constraints.  The  number  of  CFD  calls  is 
used  as  the  figure  of  merit  since  the  CFD  solutions  are  by  far  the  costliest  part  of  the 
analysis,  and  since  this  measure  is  independent  of  the  type  of  computer  used.  The 
CFD  calls  for  gradient  calculations  are  excluded  because  there  are  methods  of  finding 
the  gradients  analytically  or  semi-analytically,  without  relying  on  finite-differencing 
the  results  of  CFD  solutions,  as  was  done  in  this  research.  (See  Ref’s.  4,  13.  and 
14.) 
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For  all  these  optimizations,  the  initial  values  of  a  and  were  18.0  deg.  and  3.74 
times  the  combustor  exit  height,  respectively.  (This  corresponds  to  a  rarnp  cun'ature 
coefficient,  of  0.005.)  Tl»e  forebody,  compression  ramp,  and  inlet  configura¬ 

tions  were  held  constant  at  the  design  shown  in  Fig.  3.4.  The  forebody  and  cowl 
drag  and  the  combustor  inlet  properties  were  held  constant.  However,  the  contribu¬ 
tion  of  drag  due  to  base  pressure  was  allowed  to  var>  as  the  nozzle  geometry  changed. 
(If  the  height  of  the  nozzle  exit  plane  wa.s  less  than  the  height  of  the  forebody,  then  the 
pressure  on  the  bluff  base  was  set  to  zero.  (See  Fig.  4.1.)  TTic  nozzle  height  was  not 
allowed  to  be  greater  than  the  forcbody  height  for  the  nozzle  optimizations,  and  the 
nozzle  length  was  fixed.  The  only  constraint  placed  on  the  solution  was  that  the  pres¬ 
sure  at  the  tip  of  the  ramp  should  be  at  least  0.41  times  the  freestream  pressure.  This 
constraint  was  chosen  to  prevent  reverse  flow  into  the  nozzle  from  the  freestream,  and 
the  limiting  value  of  0.41  was  rather  arbitrarily  chosen  by  examining  the  data  from 
Ref.  4. 

The  following  side  constraints  were  imposed  on  the  optimizations:  the  initial 
ramp  angle,  a,  was  to  be  between  10  and  26  degrees;  the  nozzle  exit  plane  height  was 
to  be  between  4.0  and  6.5815  (which  was  the  foiebody  height).  No  direct  constraints 
were  placed  on 

For  the  Taylor  series  and  the  GLA's,  move  limits  of  10%  were  arbitrarily 
imposed  on  the  design  variables.  That  is,  for  each  iteration,  the  design  variables  were 
limited  to  be  within  10%  of  the  baseline  design  for  that  iteration. 

4.4.3  Optimization  Results 

The  results  for  all  the  nozzle  optimizations  aic  shown  in  TABLE  4.2.  Shown  in 
the  table  are  the  final  values  of  a  and  net  thrust  (normalized  by  p„cl,L„p  where 

=s  1.0  meters.),  and  the  required  number  of  caMs  to  the  CFD  .solver.  For  consis¬ 
tency,  the  optimized  net  thrust  values  shown  were  all  were  evaluated  by  applying  the 
CFD  solver  to  the  optimum  design  found  by  each  mcchod.  Also  shown,  for  each 
approximation  concept,  is  the  cost  in  c./mputer  time  for  1  step  in  an  iteration,  normal¬ 
ized  by  the  cost  for  1  step  using  CFD  only.  Fig.  4.6  shows  the  nozzle  contours  for  the 
initial  and  all  optimized  designs. 

Fig.  4.6  shows  that  all  the  optimized  designs  led  to  relatively  sharp  afterbody 
trailing  edges,  in  order  to  eliminate  base  drag.  Thus  the  differences  in  the  various 
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TABLE  4.x  Risulh  of  noxzle/aft«rbady  opti»nizutions.  (All  net  thrust  values  are  as 
fount'  by  CFD  solver.) 


*  As  calculated  by  the  CFD  solver. 

optimized  designs  have  to  dc  with  how  the  optimizer  varied  the  ramp  contour  between 
the  fixed  throat  location  and  the  essentially  fixed  trailing  edge. 

TABLE  4.2  shows  that  the  CFD-only  optimization  required  22  cjills  to  the  CFD 
solver  (exclusive  of  calls  for  calculating  gradients).  Each  of  these  calls  required  66 
minutes  of  CPU  time  on  an  IBM  RS  6000  workstation. 

The  uncGirected  1-D  and  MOC  optimum  designs  were  very  different  from  the 
CFD  optimization.  They  both  optimized  to  the  same  design,  because  each  hit  the 
maximum  allowable  a,  and  then  set  the  base  area  nearly  to  zero  by  adjusting  (See 
Fig.  4,6c.)  In  addition,  their  optimum  values  of  were  6%  lower  than  for  the 
CFD,  Taylor  scries  and  GLA  methods. 

The  Taylor  series  and  the  MOC  GLA  optimized  to  design.^  very  close  to  that 
which  the  CFD-only  optimization  selected,  but  did  so  with  only  7  calls  to  the  CFD 
solver  respectively.  This  represents  a  savings  in  computer  time  of  68%.  Although 
10%  move  limits  were  imposed  on  the  design  variables  for  these  optimizations,  it  is 
possible  that  greater  move  limits  could  be  imposed,  perhaps  resulting  in  some  further 
savings.  Although  the  1-D  GLA  optintum  design  is  not  quite  as  close  to  the  CFD  opti- 
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rauin  design,  the  optimum  net  thrust  is  almost  identical,  and  the  savings  in  CFD  calls 
are  the  same  as  for  the  Taylor  series  and  MOC  GLA. 

The  relative  cost  of  one  step  within  a  given  iteration  was  about  23%  cheaper  for 
the  1-D  GLA  than  for  the  MOC  GLA.  It  was  about  15%  cheaper  for  the  I-D  GLA 
than  for  the  Taylor  scries.  Each  of  these  methods  required  about  the  same  number  of 
steps  within  each  iteration  (about  14).  Thus,  if  the  minor  error  in  the  1-D  GLA  solu¬ 
tion  is  acceptable,  it  is  the  most  cost-effective  of  the  approximation  concepts  tested  in 
the  nozzle  region. 

Optimization  histories  for  the  CFD-only,  Taylor  series  and  the  GLA’s  are  shown 
in  Fig.  4.7  through  Fig.  4.9.  The  effectiveness  of  the  Taylor  series  and  the  GLA’s  are 
highlighted  by  Fig.  4.7,  which  portrays  versus  CFD  calls.  Fig.  4.8  and  Fig.  4.9 
depict  the  histories  of  a  and  respectively.  Note  that  the  Taylor  scries  and  1-D 
GLA  curves  arc  relatively  smooth  up  to  the  6nal  design,  whereas  the  MOC  GLA  his¬ 
tories  aie  rather  jagged.  But  all  these  methods  find  nearly  the  same  optimum  design. 


4.3  Summary 


The  Taylor  seric.’  local  approximation  and  the  1-D  GLA  and  MOC  GLA  were 
seen  to  provide  suitable  substitutes  for  the  CFD  solver  in  optimizing  the  nozzle.  They 
all  found  nearly  tlie  same  optimum  design  as  the  CFD-only  optimization  did,  but  v/irh 
a  savings  in  comptiter  resources  of  68%.  Their  performance  was  presaged  by  the 
single-variable  paraincterizations,  in  which  they  all  followed  the  trends  of  the  CFD 
solver,  even  when  far  from  the  design  point  upon  which  the  approximations  were 
based.  Note  that  the  neither  the  unconccted  1-D  nor  the  uncorrected  MOC  finds  the 
same  optimum  design  as  the  CFD  does. 

Thus,  in  regions  in  which  the  Sow  is  basically  isentropic,  the  approximation  con¬ 
cepts,  as  currently  used  in  structural  optimization,  only  require  modest  modifications 
to  be  applied  to  this  aerodynamic  optimization  problem.  The  savings  in  computer 
resources  would  be  much  greater  if  a  larger  number  of  design  variables  were  used. 


f 
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(a)  Initial  design  contour. 
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(b)  Optimum  design  found  by  CFD,  Taylor  and  GLA’s. 
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(c)  Optimum  design  found  by  uncorrected  1-D  and  MOC. 


Fig.  4.6  Nozzle  contours  for  the  initial  coneguration 
and  for  the  optin^um  designs  found  by  various  methods. 


Fig.  4.7  Optimization  history  of  for  CFD-only,  Taylor 
series,  and  GLA’s. 


Fig.  4.8  Optimization  history  of  a  for  CFD-only,  Taylor 
series,  and  GLA’s. 
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Fig.  4.9  Optimization  history  of  for  CFD-only, 
Taylor  series,  and  GLA’s. 


Chapter  5 


Foiebody/lnlet  Optimizations 


5. 1  Introduction 

Having  found  the  approximation  concepts  to  be  successful  in  the  no/zle  region, 
where  the  flow  near  the  walls  is  essentially  shock-free,  these  techniques  were  next 
applied  to  the  optimization  of  the  forebody/inlet  region.  The  flow  there  is  very  differ¬ 
ent  from  the  nozzle  flow,  as  it  is  characterized  by  the  presence  of  oblique  shock 
waves.  It  is  not  clear  beforehand  how  there  shock  waves  might  affect  the  approxima¬ 
tion  concept  techniques.  Thus  there  is  a  great  deal  more  risk  involved  in  trying  to 
apply  these  techniques  there. 

This  chapter  describes  the  formulation  of  the  design  problem,  in  which  the  opti¬ 
mizer  is  allowed  to  choose  designs  without  explicit  regard  to  the  resulting  shock 
structures.  (This  is  in  contrast  to  actual  design  practice,  which  only  uses  “properly” 
placed  shocks.)  It  then  describes  the  special  considerations  required  because  of  the 
presence  of  imperfectly-placed  shocks.  Next  it  outlines  the  selection  of  suitable 
design  variables.  Then  the  results  of  the  optimizations  are  presented. 

5.2  General  Design  Considerations 

The  2-D  scramjet  vehicle  uses  a  mixed-compression  inlet.  Such  an  inlet  can  be 
constructed  by  combining  a  nose/  forebody,  compression  ramp  and  cowl,  as  shown  in 
Fig.  5.1.  Ifhe  nose/forebody  region  is  a  straight  wedge,  whose  oblique  shock  pro¬ 
vides  the  initial  compression  for  the  scramjet  engine.  The  compression  ramp  gener¬ 
ates  a  second  shock  to  provide  further  compression.  Ideally,  these  two  oblique  shocks 
would  impinge  on  the  cowl  lip.  A  third  (internal-compression)  shock  is  formed  at  the 
cowl  lip  and  ideally  would  impinge  upon  the  corner  of  the  compression  ramp,  as 
shown  in  Fig.  5.1.  Thus  shock  would  then  interact  with  and,  ideally,  would  exactly 
cancel  the  expansion  fan  created  at  the  ramp  comer.  In  such  an  ideal,  design,  the  com¬ 
bustor  is  provided  with  uniform, '1-D  flow  at  high  static  pressure.  Ifhe  pressure  rises 


Combustor 
Inlet  Plane 


Fig.  5.1  Forebody/Inlet  geometry. 


across  the  nose  and  compression  ramp  shocks  also  create  a  ‘'forebody  drag”.  The 
cowl  shock  can  also  contribute  to  fcrebody  drag  if  it  is  not  properly  placed  at  the  ramp 
comer.  There  is  also  a  small  contribution  to  drag  from  the  oblique  shock  on  the  lower 
surface  of  the  cowl  lip. 

For  ideally-placed  shocks  and  given  freestream  conditions,  the  design  problem 
then  becomes  one  of  determining  the  nose  and  compression  ramp  angles,  and 
which  strike  the  best  balance  between  high  combustor  performance  and  low 
forebody  drag.  For  given  and  the  parameters  which  determine  whether 
the  shocks  arc  ideally-placed  are:  the  forcbody  length,  the  compression  ramp 
length,  the  inlet  cowl  length,  and  the  inlet  height 

In  this  research,  to  demonstrate  the  use  of  approximation  techniques  for  use  in 
design  optimization,  the  following  problem  was  chosen.  Given  an  initial  design  with 
non-ideally  placed  shocks  and  arbitrary  0,,,,  and  0„„p,  optimize  to  the  design  which 
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maximizes  f  Thus,  v/e  expect  the  optimizer  to  find  the  best  and  and  to 
place  the  shocks  to  best  advantage.  (As  explained  in  Section  5.3,  this  will  not  neces¬ 
sarily  correspond  to  perfectly-placed  shocks,  if  the  effects  of  flow  nonuniformity  in 
the  inlet  duct  are  not  taken  into  account.) 

The  general  utility  of  this  approach  is  that,  for  a  more  complicated  and  realistic 
design  problem  (such  as  for  a  3-D  vehicle  with  viscous,  reacting  flows)  it  may  be  dif¬ 
ficult  to  specily,  a  priori,  a  configuration  which  will  perfectly  place  the  shocks.  Thus, 
a  method  which  would,  as  a  natural  consequence  of  the  optimization  process,  place 
the  shocks  in  the  best  way,  would  be  a  powerful  design  tool.  The  approximation  tech¬ 
niques,  such  as  Taylor  series  and  the  Global-Local  Approximations  may  hold  this 
promise,  and  the  problem  as  stated  in  the  previous  paragraph  is  one  way  of  assessing 
their  validity.  More  specifically,  if  an  approximation  technique  can  be  made  to  opti¬ 
mize  to  the  same  configuration  as  the  CFD,  then  that  technique  would  become  an 
attractive  alternative  to  using  the  CFD  alone  in  the  optimization  process. 

The  following  sections  describe  how  the  various  modeling  levels  and  approxima¬ 
tion  techniques  were  applied  to  this  problem. 


5.3  Special  Design  Considerations 

The  presence  of  shocks  which  impinge  upon  the  walls,  and  which  may  interact 
with  each  other  and  with  expansions  in  the  flow  makes  this  problem  essentially  very 
different  from  the  problem  of  optimizing  the  nozzle.  Some  discussion  of  the  special 
problems  encountered  in  this  design  problem  will  serve  to  make  clear  the  approach 
that  was  used  to  perform  the  optimizations,  and  their  results. 

In  practice,  it  is  undisputed  that  Che  only  desirable  shock  configuration  for  the 
mixed-compression  inlet  is  perfect  shock  placement,  as  described  in  the  introduction 
to  this  chapter.  However,  given  that  this  research  is  limited  to  inviscid  flows,  without 
consideration  of  shock/boundary  layer  interactions,  and  without  any  details  of  the 
flow  in  the  combustor,  some  of  the  impetus  to  place  the  shocks  perfectly  is  missing. 
Thus,  it  is  possible  that  the  optimizer  will  select  non-perfect  shock  placements.  A  dis¬ 
cussion  of  why  this  might  be  so  is  warranted. 
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5.3.1  Effects  of  Non>Perfect  Shock  Placement 

Neglecting  viscous  and  unsteady  effects,  the  consequences  of  imperfect  placing 
of  thb  shocks  ju*e  as  follows. 

a. )  If  the  cowl  shock  hits  the  compression  ramp  short  of  the  corner,  as  shown  in 
Fig.  5.2a,  a  reflection  results.  (This  shock  may  be  strong  eriough  to  become  a  normal 
shock,  thus  unstarting  the  engine.)  The  very  high  pressure  jump  downstream  of  this 
reflection  causes  an  increase  in  forebody  drag,  which  can  be  quite  severe.  In  addition, 
the  reflected  shock  interacts  with,  but  does  not  cancel  the  expansion  fan  from  the  cor¬ 
ner.  The  result  is  a  complicated  system  of  shocks  and  expansions  in  the  inlet  duct, 
with  a  resulting  flow  nonuniformity.  This  non-uniforraity  can  be  extremely  detrimen¬ 
tal  to  the  combustion  process.  The  only  positive  aspect  of  having  the  cowl  shock  fall 
short  of  the  ramp  comer  (assuming  it  does  not  unstart  the  engine)  is  that  it  provides  an 
additional  compression,  which  may  contribute  to  improved  combustor  performance,  if 
the  nonuniformity  is  not  too  great, 

b. )  If  the  cowl  shock  passes  inside  the  inlet  duct,  as  shown  in  Fig.  5.2b,  there  is 
again  an  non-cancelling  interaction  with  the  expansion  fan,  and  a  resulting  flow  non- 
uniformity  in  the  combustor  entrance.  However,  there  is  no  additional  forebody  drag. 

c. )  If  the  nose  shock  falls  short  of  the  cowl  lip,  as  shown  in  Fig.  5.2c,  the  nose 
section  of  the  forebody  is  longer  than  it  needs  to  be,  with  a  consequent  increase  in 
forebody  drag,  and  no  more  engine  mass  flow  than  if  the  shock  were  perfectly  placed. 
In  addition,  the  flow  around  the  cowl  lip  is  changed,  since  the  angle  of  attack  there  is 
now  nonzero.  This  will  result  in  a  reduction  in  cowl  drag,  since  the  effective  wedge 
angle  on  the  upper  surface  is  decreased.  Due  to  the  small  vertical  area  of  the  cowl, 
diis  reduction  in  cowl  drag  does  not  usually  offset  the  forebody  drag  increase.  So  the 
net  result  of  having  the  nose  shock  fall  short  of  the  cowl  is  usually  a  reduction  in  the 
net  thrust.  Similar  changes  occur  if  only  the  compression  ramp  shock  or  both  the  nose 
and  compression  ramp  shocks  fall  short  of  the  nose.  Thus,  there  will  be  a  natural  ten¬ 
dency  for  the  optimizer  to  avoid  solutions  which  place  these  shocks  short  of  the  lip. 

d. )  If  the  nose  shock,  compression  ramp  shock,  or  both  are  too  shallow,  and  hit 
inside  the  cowl  (see  Fig.  5.2d),  they  would  interact  with  and  cause  changes  in  the 
strength  and  angle  of  the  cowl  shock.  This,  in  turn,  could  cause  improper  placement 
of  the  cowl  shock,  with  all  of  the  problems  shown  in  a.)  and  b.)  above.  But,  it  is  pos¬ 
sible  to  have  the  nose  and/or  ramp  shock  hit  inside  the  cowl,  and  have  the  resulting 
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(a)  Cowl  shock  hits  short  and  reflects  off  ramp. 


(b)  Cowl  shock  enters  inlet  duct  directly. 


(c)  Nose  shock  misses  cowl. 


(d)  Nose  and  cowl  shocks  hit  inside  cowl. 
Fig,  5.2  Improperly-placed  inlet  shocks. 


cowl  shock  be  perfectly-placed,  or  nearly  so.  This  indeed  happened  on  some  of  the 
optimizations  to  be  presented  in  Section  5.5.  In  fact,  it  was  observed  that,  during  the 
optimization  process,  once  a  shock  has  been  placed  so  as  to  impinge  inside  the  duct,  it 
is  almost  impossible  to  get  it  back  outside  again. 

53.2  Inlet  Duct  Flow  Non<Uniforniity 

The  CFE)  solver  naturally  accounts  for  all  of  the  shock-shock  and  shock-expan¬ 
sion  interactions  in  the  inlet  duct.  It  then  averages  the  flow  properties  in  the  combus¬ 
tor  inlet  plane,  and  passes  this  information  to  the  1-D  combustor  model.  On  the  other 
hand,  the  approximate  flow  model  for  the  forebody/inlet  region  terminates  at  the  cor¬ 
ner  of  the  compression  ramp  (i.e.  in  the  inlet  plane),  and  averages  the  properties  there. 
(See  Section  3.4.)  Thus  this  model  does  not  account  for  the  shock-expansion  interac¬ 
tions  and  reflections  inside  the  duct.  These  effects  influence  both  the  flow  non-uni¬ 
formity  at  tlie  combustor  inlet,  as  well  as  the  averaged  properties  being  passed  to  the 
combustor  model.  This  is  a  major  reason  that  the  CFD  solver  is  herein  considered  to 
be  more  accurate  than  the  approximate  model. 

In  this  research,  no  penalty  was  assessed  to  the  combustor  performance  for  such 
flow  non-uniformity.  It  is  not  a  simple  matter  to  determine  such  a  penalty.  In  addi¬ 
tion,  as  the  primary  purpose  of  this  work  is  to  demonstrate  approximation  techniques 
which  accurately  mimic  the  CFD  solutions  in  the  optimization  process,  this  addi¬ 
tional  complexity  is  not  required,  as  long  as  the  optimum  designs  do  not  posses  unrea¬ 
sonably  high  non-uniformities. 

An  alternate  procedure  for  handling  flow  non-uniformity  is  to  impose  a  con¬ 
straint  in  the  optimization  procedure.  This  was  tried  by  using  the  approximate  model 
and  optimizing,  imposing  the  constraint  that  the  standard  deviation  of  the  pressures 
(at  the  CFD  grid  points)  in  the  inlet  plane  be  less  than  a  certain  value.  Depending 
upon  which  set  of  design  variables  was  used  (Sec  Section  5.4),  the  optimizer  did  one 
of  the  following:  either  froze  and  (if  not  directly,  then  by  whatever 

design  variables  control  them)  and  adjusted  the  remaining  design  variables  to  place 
the  cowl  shock  perfectly  at  the  ramp  comer,  then  stopped;  or  reduced  and 
as  low  as  possible  and  then  sought  an  optimum  there.  Both  of  these  results  are  reflec¬ 
tions  of  the  fact  that  the  optimizer  places  first  prit»rity  on  insuring  that  the  design  is 
fea;>ible  (i.e.  not  in  violation  of  any  constraints).  It  will  do  so  to  the  exclusion  of 
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improving  the  objective  function,  if  necessary,  until  the  design  is  feasible.  And  since 
the  amount  of  non-uniformity  increases  very  rapidly  as  we  move  away  from  perfect 
cowl  shock  placement,  the  optimizer  finds  a  very  steep  local  optimum,  from  which  it 
is  unlikely  to  extract  itself. 

For  these  reasons,  no  constraints  on  inlet  duct  flow  non-uniformity  were  imposed 
in  this  research,  nor  were  any  penalties  imposed  on  the  combustor  performance  for 
having  such  non-uniformity.  However,  this  is  somewhat  mitigated  by  the  fact  that 
certain  influences  will  tend  to  keep  the  cowl  shock  from  being  grossly  misplaced. 
First,  if  this  shock  hits  short  of  the  ramp  comer  (  Fig.  5.2a),  the  reflection  causes  an 
increase  in  forebody  drag,  which  can  offset  any  increase  in  thrust.  Secondly  if  the 
shock  passes  inside  the  duct,  as  in  Fig.  5.2b,  the  average  pressure  in  the  inlet  plane  is 
lowered,  reducing  combustor  performance.  Together,  these  two  influences  act  to 
keep  any  misplacement  of  the  cowl  .shock  (and  the  resulting  flow  non-uniformities) 
from  becoming  wildly  excessive. 


5  J.3  Suiiiinai7  of  Design  Considerations 

No  satisfactory  way  was  found  to  explicitly  prevent  the  imperfect  placing  of  the 
shocks,  or  of  penalizing  the  combustor  performance  for  flow  non-uniiformity  in  the 
inlet  duct.  However,  certain  influences  were  seen  to  dictate  against  gross  non-unifor¬ 
mities  and  against  having  the  nose  and  cowl  shocks  fall  short  of  the  cowl  lip.  There 
is  no  such  in^ucnce  to  keep  the  nose  and/or  ramp  shock  from  entering  the  duct,  as 
long  as  this  does  not  cause  the  cowl  shock  to  be  grossly  misplaced. 

5.4  Choice  of  Design  Variables 

As  in  the  nozzle  optimizations  of  Chapter  4,  an  investigation  was  made  to  find 
the  best  set  of  design  variables  to  use  in  optimizing  the  forebody/inlet.  To  do  this,  the 
approximate  model  (based  on  oblique  shock  theory  for  the  forebody/inlet,  and  the 
Method  of  Characteristics,  or  MOC,  for  the  nozzle)  was  used  in  the  optimizer,  and 
optimizations  were  run  for  various  sets  of  design  variables.  The  resulting  optimizer 
behavior  was  ob.scrved.  The  design  variable  .set  which  led  to  the  best  behavior  was 
then  selected  as  the  best  set. 
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The  followiiig  design  variable  sets  were  considered.  Their  meanings  are 
explained  in  Fig.  5.3. 

^ramp’  ^/hi  ^ramp'  ^tld 
^-)  ^noif>  '^iioj«»  ^fb'  ^ramp*  ^cawt 

^noje»  ^fiote'  ^rAmp'  ^ramp'  ®tld  I'cowt 
tl<)  ^note^  ^^ramp'  ^^ramp*  ®*td  ^>co*l 

Design  variable  set  a.)  is  the  most  intuitive  one,  because  it  is  these  variables  which 
directly  control  the  shock  strengths  and  angles,  and  the  placement  of  the  shocks. 
However,  there  is  a  disadvantages  to  using  this  set  in  the  optimization,  Let  us  exam¬ 
ine  the  behavior  of  the  gradient  calculations  (for  finding  the  search  direction)  using 
set  a.).  Assume  an  initial  design  that  has  nearly  f>erfect  shock  placement.  We  calcu¬ 
late  the  gradient  clement  b>  incrementally  increasing  0^^^^  while  leaving  all 

FO.!Kp 

Other  design  variables  constant.  In  so  doing,  we  might  expect  to  see  an  increase  in 
combustor  performance  due  to  the  increased  ramp  compression.  We  would  also 
expect  an  increase  in  forebody  drag,  both  due  to  the  increase  in  pressure  on  the  ramp 
(which  is  unavoidable),  but  also  due  to  the  increase  in  the  ramp  height.  (See  Fig. 
5.4a.)  In  addition,  since  the  origin  of  the  ramp  shock  has  not  changed,  the  increase  in 
oblique  shock  angle  may  cause  the  ramp  shock  to  miss  the  cowl.  The  increased  ramp 
height  and  improperly-placed  shock  may  nullify  the  improvement  in  ramp  compres¬ 
sion.  Accordingly,  the  gradient  calculation  would  indicate  that  there  is  no  benefit  to 
be  had  by  increasing  9^^^. 

The  natural  solution  to  this  is  to  decrease  as  we  increase  This  would 

have  the  effects  of  reducing  the  ramp  height,  and  moving  the  ramp  shock  closer  to  the 
cowl  lip.  Thus,  we  might  increase  inlet  compression  with  a  reduced  increase  in  fore¬ 
body  drag,  while  simultaneously  helping  to  maintain  a  favorable  shock  placement. 
We  might  then  find  that  there  is  actually  a  great  deal  to  be  gained  by  increasing 
However,  this  logic  is  missing  in  a  gradient  calculation  performed  with  design  vari¬ 
able  set  a.).  Thus,  a  beneficial  search  direction  may  well  be  not  be  recognized  by  the 
optimizer.  (Indeed,  in  the  tests  of  the  design  variable  sets,  set  a.)  often  caused  the  opti¬ 
mizer  to  fail  to  see  beneficial  search  directions  which  would  be  obvious  to  an  engi¬ 
neer.) 


60 


So,  it  seems  sensible  that  we  should  build  into  the  gradient  calculations  the  fact 

that,  in  reality,  we  would  often  change  some  of  the  design  parameters  simultaneously. 

Design  variable  sets  b.)  through  d.)  were  selected  because  they  build  ^uch  multi- 

vaiiable  flavor  into  the  gradient  calculations.  Thus,  they  may  be  considered  “smart.  ” 

dt  * 

design  variable  sets.  Consider,  for  example,  .set  b.).  If  we  calculate  ,  we  will 
incrementally  increase  h^,,  while  holding  (as  well  as  the  nose  design* variables) 
constant,  which  effectively  requires  shortening  (Sec  Fig.  5.4b.)  For  a  given 

increase  in  the  value  of  L^„„p  will  be  less  than  it  would  be  when  using  design 
variable  set  a.).  This  also  implies  that  is  shorter.  And  so  the  optimizer  will,  in 
one  element  of  the  giadient  calculation,  consider  a  search  direction  which  simulta¬ 
neously  seeks  the  benefits  of  increased  ramp  compression  and  decreased  forebody 
drag,  while  making  an  effort  to  maintain  the  proper  shock  placement. 

TABLE  5.1  shows  the  results  of  the  design  variable  tests.  Note  that  the  initial  condi¬ 
tions  for  all  o.f  tliesc  Tests  (when  converted  into  terms  of 

are  as  foiiow's;  8.0  deg,  18.0  deg,  =  35.0,  =  13.0,  and 

"  2.0,  .Note  also  that  the  nozzle  design  was  the  same  as  the  initial  condition  used 
in  the  nozzle  tests  of  Chapter  4,  with  a  =  18.0  deg,  and  =  0.005.  The  inlet  height 
was  fixed  at  0.6667.  Tiie  free  stream  flow  conditions  were  those  at  =  6.0  and  30 
km  altitude.  The  maximum  allowable  ramp  angle  was  20  0  deg. 

In  T.4BLE  5. 1  the  optimum  design  for  each  of  the  design  variahln  sets  is  shown, 
as  well  as  the  resulting  net  thr  ust.  (Note  that  these  net  thrust  values  are  those  found 
by  the  approximate  model,,  and  not  by  the  CFD  solver.)  From  these  tests,  it  was  clear 
that  the  design  variable  set  b.),  that  is  and  is  the  most 

attractive  one.  Set  a.)  did  not  allow  the  ramp  angle  to  increase,  for  the  reasons  men¬ 
tioned  previously  in  this  section.  Set  c.)  made  the  nose  and  ramp  angles  decrea.se, 
counter  to  reason.  Although  set  d.)  found  reasonable  nose  and  ramp  angles,  the  cowl 
shock  was  badly  placed,  passing  directly  inside  the  inlet  duct,  instead  of  impinging 
on  the  ramp.  (Similar  to  the  ramp  shock  configuration  shown  in  Fig.  5.2b.)  Based 
on  its  superior  performance  in  these  tests,  design  variable  set  b.)  was  used  for  all  sub-^ 
sequent  optimizations  of  the  forebody/inlet  region. 


ParturtMKJ 

Oe$i(|n 


(a)  Using  design  variable  set  a.). 


Design 


(b)  Using  design  variable  set  b.)- 


Rg,  5,4  EfiFects  on  search  direction  calculation  of 
choosing  difiisrent  design  variable  sets. 
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TABLE  S.1 _ 

Results  of  design  variable  tests:  optimum  configurations  for  each  set  of  design 
variables.  (All  design  variables  are  converted  into  angles  and  lengths.) 


5.5  Optimization  Results 

Optimizations  of  the  forebody/inlet  design  were  performed  using  each  of  the  fol¬ 
lowing  modeling  levels  and  approximation  techniques:  CFD  only;  approximate 
model  (based  on  oblique  shock  theory  for  the  forcbody/inlet,  and  MOC  for  the  noz¬ 
zle);  Taylor  series;  GLA  based  on  the  approximate  model.  Just  as  for  the  nozzle  opti¬ 
mizations,  the  objective  function  was  the  normalized  net  thrust,  which  has  been 
normalized  by  Side  constraints  were  imposed  to  keep  the  design  variables 

within  reasonable  limits,  and  inequality  constraints  were  placed  on  and  0^^^^  to 
keep  them  within  certain  limits  (6.0^0,^^,  s  11.0  and  11.5  S0^^^s  20.0).  The  inlet  height 
was  fixed  at  0.6667. 

The  initial  conditions  for  these  optimizations  were  (converted  from  the  actual 
design  variables  into  those  which  are  more  easily  visualized): 

Ko.,  =  8-0 

=  18'0deg 
=  35.0 
i-™*,  =  13.0 
=3.20 

To  gain  insight  into  what  the  optimum  solutions  should  generally  look  like,  a  pre¬ 
liminary  optimization  was  run  with  a  special  program  which  always  placed  the  shocks 
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perfectly,  and  which  used  only  the  approximate  model.  Since  the  shock  placcraeni  in 
such  cases  is  uniquely  determined  by  and  these  were  the  only  two  design 
variables.  The  optimum  design  from  this  test  selected  the  maximum  value  of  9^^^^ 
(20.0  deg),  and  very  nearly  the  maximum  value  of  9^„,  (10.9  deg).  This  is  because 
the  inlet  compression  has  not  yet  reached  that  which  is  most  efficient  for  the  combus¬ 
tor.  (See  Section  3.3  and  Fig.  3.9.)  It  is  likely  that  9,^^^,  would  have  been  increased 
even  further,  had  the  inequality  constraints  not  prevented  that. 

Thus,  we  should  expect  that  the  optimizations  to  be  performed  with  non-ideal 
shock  placement  should  generally  follow  the  trend  of  seeking  increased  ccrapression 
from  both  the  nose  and  ramp  angles. 

5.5.1  Performance  of  the  GLA 

It  was  observed  early  on  that  the  GLA  method  behaved  very  erratically  within 
any  given  iteration.  Fig.  5.5  shows  the  history  of  (on  a  logarithmic  scale)  for 
the  steps  within  the  second  iteration.  (Recall  that  an  iteration  is  a  complete  optimiza¬ 
tion  performed  using  only  the  corrected  approximate  method,  and  based  on  a  given 
baseline  design.  A  step  is  a  single  function  evaluation  within  that  optimization.)  The 
uncorrected  approximate  method  history  is  shown  for  comparison,  since  it  was  found 
to  be  a  reasonably  good  approximation  to  the  CFD  solver.  (See  TABLE  5.2)  The 
figure  shows  the  very  erratic  behavior  of  the  GLA.  The  values  of  shown  are  as 
found  by  the  approximate  method  and  by  the  GLA  respectively,  and  have  not  been 
evaluated  by  the  CFD  .solver. 

The  wild  variations  in  the  GLA  history  were  found  to  be  due  in  large  part  to  diffi¬ 
culties  that  this  method  has  in  dealing  with  shifting  shock  impingement  points  on  the 
ramp  surface  near  the  inlet.  Basically,  as  a  shock  moves  to  a  new  position,  the  correc¬ 
tion  factors  of  the  GLA  do  not  account  for  this,  and  improperly  magnify  or  reduce  the 
pressures  in  the  vicinity  of  the  shock.  In  addition,  the  fact  that  the  CFD  solver  and 
the  approximate  method  do  not  capture  shocks  equally  crisply  causes  further  en  or. 

The  GLA  and  Taylor  series  methods,  as  heretofore  implemented,  share  an  inabil¬ 
ity  to  account  for  the  effect  of  discrete  phenomena  such  as  shocks  that  shift  position  as 
design  variables  are  changed.  This  effect  is  most  crucial  when  it  is  manifested  in  the 
calculation  of  correction  factors  for  the  GLA,  and  somewhat  less  critical,  although 
still  important,  if  it  occurs  in  the  calculation  of  the  search  direction,  for  either  method. 
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Note;  10.0  is  added  to  to  allow  log  plotting  of  neg  values. 


Fig.  5.5  Behavior  of  GLA  and  approximate 
methods  within  one  iteration.  (Note  log  scale.) 

To  illustrate  this  effect,  consider  one  element  of  the  wall  pressure  gradient  calcu¬ 
lation  for  a  GLA  which  will  involve  a  shock  which  moves  during  the  calculation. 
(See  Fig.  5.6.)  Let  Xj  be  the  design  variable  of  interest.  Fig.  5.6a  shows  the  shock 
impingement  point  moves  from  point  A  to  point  B  as  the  design  variable  is  perturbed 
for  the  derivative  calculation.  Fig.  5.6b  shows  the  wall  pressure  distributions  for  the 
two  shock  configurations.  Note  that  the  change  in  pressure  is  much  greater  in  region 
n  than  in  either  of  regions  I  or  in.  This  is  manifested  in  the  derivative  calculation 
plotted  in  Fig.  5.6c,  in  which,  for  compactness  in  the  figure,  the  peak  value  is  actu¬ 
ally  very  much  under-represented,  since  the  denominator  of  the  calculation  is  very 
small,  due  to  the  small  perturbation  in  the  design  variable.  Thus,  for  calculation  of 
correction  factors,  or  for  calculating  a  search  direction,  the  gradient  calculation  dic¬ 
tates  that  the  effect  of  the  shock  movement  is  to  be  dominated  by  an  extreme  sensitiv¬ 
ity  to  the  design  variable  in  region  II  (which  is  fixed  in  space). 

In  reality,  we  know  that  the  spike  in  region  II  is  an  artifact  of  not  recognizing  that 
the  shock  is  in  transit.  A  more  reasonable  representation  of  the  derivative  is  that  in 
Fig.  5.6d,  in  which  there  is  a  moving  shock  location,  and  separate  derivatives 
upstream  and  downstream  of  the  shock. 


B  A 

(l)  Shock  locations  for  initial  and  perturbed  design  variable. 


c.)  Derivative  of  pressure  wrt  design  variable. 


Pre-shock 


Postshock 


d.)  Actual  derivative,  if  properly  calculated. 


Fig.  5.6  Effect  of  moving  shock  impingement  point  in 
gradient  calculation. 
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Interestingly,  the  inability  to  account  for  the  movement  of  the  shock  is  somewhat 
alleviated  when  using  CFD  calculations  for  performing  the  derivatives.  This  is 
because  the  shock  dissipation  inherent  in  the  CFD  solution  “spreads  out”  the  shock, 
thus  giving  adjacent  points  some  warning  that  the  shock  is  approaching,  This  par¬ 
tially  explains  why  the  Taylor  series  approach  is  better  behaved,  if  not  more  accurate, 
than  the  GLA.  In  addition,  the  Taylor  series  involves  only  one  set  of  gradient  calcu¬ 
lations,  namely  those  based  on  the  CFD.  The  GLA  uses  gradients  based  both  on  the 
CFD  and  on  the  approximate  method,  which  uses  sharp  shock  discontinuities,  such  as 
those  in  Fig.  5.6.  The  problem  is  thus  magnified  for  the  GLA. 

The  obvious  solution  to  this  problem  is  to  calculate  separate  gradient  and  correc¬ 
tion  factor  calculations  for  pre-  and  post-shock  locations,  and  apply  them  only  where 
appropriate.  This  requires  devising  two  adjoining  non-dimensionalized  coordinate 
systems  for  each  wall  surface,  with  the  physical  location  of  the  interface  to  be  deter¬ 
mined  (by  the  approximate  method)  once  the  location  of  the  shock  is  known.  This 
process  is  depicted  in  Fig.  5.7.  As  the  optimization  is  performed  using  the  GLA,  at 
each  step  the  shock  location  is  first  found  using  the  approximate  method,  and  the 
approximate  method  wall  pressures  are  then  corrected  by  application  of  the  pre-  or 
post-shock  correction  factors,  as  appropriate.  Similar  calculations  must  be  performed 
for  the  inlet  plane  flow  properties,  since  the  intersection  of  the  cowl  shock  (or  its 
reflection)  with  the  inlet  plane  also  moves. 

One  limitation  of  this  method  is  that  it  requires  the  approximate  method  (which  is 
the  only  one  used  in  the  actual  optimizations)  to  be  able  to  predict  fairly  accurately 
the  shock  locations.  In  the  problem  at  hand,  this  is  not  a  problem,  since  the  oblique 
shock  theory  is  very  accurate  in  this  regard.  For  more  complicated  problems,  this  may 
require  more  attention. 

The  pre-and  post-shock  coordinate  system  technique  for  calculating  gradients 
and  correcting  wall  pressures  was  applied  to  one  iteration  (40  steps)  of  an  optimiza¬ 
tion  using  the  GLA.  (Note  that  the  inlet  plane  flow  properties  were  not  treated.)  The 
result  was  a  marked  decrease  in  the  erratic  behavior.  See  Fig.  5.8.  (Note  that  10.0 
has  l^een  added  to  the  values  of  .in  the  plot,  to  allow  the  negative  values  to  be 
plotted  on  the  logarithmic  scale.)  The  maximum  value  of  f„,  recorded  in  this  itera¬ 
tion  was  14.99,  compared  to  2:551  for  the  same  iteration  using  the  original  formulation 
of  the  GLA.  The  uncorrccted  approximate  model  (which  mimics  the  CFD  very  well) 
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b.)  Wall  pressures. 


c.)  Wall  pressures,  in  non>diiiiensionailzed  pre-  and  post-shock 
coordinate  systems. 


dX, 

pie-shock 

post-shock 

0  10  1 

d.)  DerivaliTe  caladations  in  non^iimensionatized  pre-  and 
portHriuKk  oftordhiate  sj^dcais. 


Bg.  S.7  Gradient  calculation  using  pre-  and  post¬ 
shock  coordinate  system. 
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Fig,  5,8  EflEcct  of  using  prc-  and  post-shock 
coordinate  systems  on  one  iteration  of  the  GLA  in 
the  forebody/inlet  region.  (Note  log  scale.) 


Fig.  5,9  Hisi(»7  of  nose  angle  for  one  GLA 
iteration,  using  prc-and  pmt-shock  cocmltnate 
systems,  with  crnnparison  to  approximate  rootkl. 
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calculated  to  be  14.96  for  the  optimum  design  in  this  iteration  of  the  new  GLA 
test.  However,  the  values  of  thrust  and  drag  found  by  the  improved  GLA  were  too 
high  (97.4  and  82.4  respectively,  versus  55.3  and  40,3  as  found  by  the  approximate 
method).  So  it  is  apparent  that  the  new  GLA  formulation  does  not  entirely  solve  the 
problem.  The  error  in  thrust  may  be  a  result  of  the  lack  of  a  pre-  and  post-shock  sys¬ 
tem  for  the  inlet  plane  flow  properties.  The  drag  error  may  also  be  due  in  part  to  the 
shock  dissipation  in  the  CFD  solution.  Because  of  this,  the  approximate  method  and 
the  CFD  do  not  capture  shocks  equally  crisply,  and  the  result  will  be  that  erroneous 
correction  factors  will  be  calculated  around  shocks.  As  the  wall  pressures  increase 
due  to  incpiasing  shock  strength  during  the  optimization,  the  errors  induced  by  the 
erroneous  correction  factors  will  grow.  Indeed,  this  increasing  error  was  found  in  the 
test  of  the  pre-  and  post-shock  coordinate  systems. 

A  history  of  the  nose  angle  for  this  iteration  is  shown  in  Fig.  5.9.  Note  that, 
although  the  erratic  behavior  of  f,,,  has  largely  been  cured  by  using  the  pre-  and 
post-shock  coordinate  systems,  the  optimizer  treats  this  model  almost  exactly  the 
same  as  it  does  the  original  GLA.  Thus,  the  poor  selection  of  an  optimum  design  by 
the  GLA  is  pertiaps  more  due  to  the  other  sources  of  error  mentioned  in  the  previous 
paragraph  than  due  to  the  need  for  a  pre-and  post-.shock  coordinate  system  Neverthe¬ 
less,  the  reduction  in  erratic  behavior  afforded  by  the  pre-  and  post  shock  system  is  an 
important  ingredient  in  formulating  a  generally-applicabie  GLA.  This  is  espccialiv 
true  when  active  behavior  constraints  (such  as  limits  on  the  pressure  at  the  ramp  up> 
are  encountered.  In  such  cases,  the  wild  fluctuations  in  wall  pressure  are  impK>rtant  in 
their  own  right 

Thus,  it  is  apparent  that  a  straightforward  application  of  GLA’s  to  the  fi>fetH»d>/ 
inlet  region,  using  the  techniqucji  used  in  structural  optimization  problems,  will  not 
ro«tt  with  success.  More  work  is  needed  in  this  area  to  unlock  whatever  potential  dvr 
GLA  method  has  for  the  forebody/inlet  region.  In  fact,  for  the  GLA  to  be  applied  to 
point  properties  in  any  flow  regime  which  has  moving,  discontinuous  phenomena  (such 
as  shocks,  transition  to  turbulence  or  .separation),  these  problems  will  have  to  be  over¬ 
come. 

Because  of  this  behavior,  a  complete  GLA  optimization,  based  on  point  proper¬ 
ties,  was  not  performed.  However,  because  the  Taylor  .scries  did  not  suffer  as  badly 
in  this  regard,  it  was  included  in  the  set  of  complete  optimizations. 
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5.5.2  GLA’s  Applied  to  Integral  Measures 

An  attempt  to  circumveni  the  problems  with  applying  the  GLA  to  point  proper¬ 
ties  was  made,  by  applying  it  to  the  integrated  properties  (thrust,  drag,  and  The 
central  idea  here  is  that  the  smoothing  process  inherent  in  integrating  would  alleviate 
the  problems  with  shock  motion. 

Several  optimizations  were  run.  It  was  found  that  there  was  no  apparent  advan¬ 
tage  to  applying  the  GLA  to  thrust  and  drag  separately,  as  opposed  to  applying  it 
directly  to  It  was  also  found  that  including  derivative  information  in  the  GLA 
caused  the  optimizer  to  select  a  different  optimum  than  the  CFD  and  approximate 
methods  did.  However,  if  the  GLA  was  limited  to  zeroth-order  (i.e.  constant  scaling 
factor),  excellent  agreement  with  the  approximate  method  (and  CFD)  was  achieved. 
It  is  important  to  note  that,  wh<;n  such  a  simple  scaling  is  applied  to  the  integral  mea- 
swes,  this  is  the  sainc  au  simply  multiplying  all  the  approximate  method  solutions 
(wtthin  one  iteration)  by  a  comsant  numter. 

SS3  OprtaintfMi  toUki 

CcMi^lcie  optiBiiaiiKNas  were  performed  using  CFD.  the  unconected  approxi- 
OMMr  mettod.  and  the  Taj^lor  series  method.  The  results  of  these  optimizations  are 
ilwwn  HI  TABLE  5.2  hkne  that  the  optimum  designs  are  presented  in  terms  of  the 
dcaign  vnnahka  of  set  a. )  w  Seebota  3.4.  even  though  the  actual  design  variables 
went  those  of  set  b.  i  This  was  done  to  simplify  intcrpreiation  of  the  resuif.. 

For  the  CFl>-only  case.  25  caib  to  the  CFD  solver  (excluding  tho.se  for  gradient 
caktdbbosM)  were  required.  The  optimized  design  increased  the  net  thrust  by  248% 
over  dm  UMtiaJ  design.  The  shock  placement  configuration  for  this  design  is  shown  in 
Fig.  5.10.  This  figvre  confirms  the  reasoning  of  Section  5.4,  in  which  it  wa.s  postu- 
laml  that  there  are  influences  which  should  keep  the  shocks  from  being  grossly  mis- 
pinevdL  Indeed,  the  improvement  in  shock  placement  is  substantial.  The  only 
significant  deviation  from  what  was  expected  based  upon  the  “perfect  shijck  place¬ 
ment'*  test  was  that  the  nose  angle  decreased  instead  of  increasing. 
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Table  S J  Results  lof  forebody^let  optimizations  (All  net  thrust  values  are  as  found  by  CFD 
solver.) 


*  As  calculated  by  the  CFD  solver. 


Note,  in  TABLE  5.2,  that  the  optimum  design  found  by  the  uncorrectrd  approximate 
model  is  close  to  that  found  by  using  only  CFD.  The  only  significant  differences  in 
design  variables  art  for  (5.6%  error)  and  (6.1%  error).  And  the  value  of 
is  8.2%  lower  than  that  found  using  only  CFD.  The  approximate  model  appar¬ 
ently  captures  all  the  relevant  physics  sufficiently  well  to  correctly  replicate  the  gen¬ 
eral  behavior  of  the  CFD  model.  And  the  cost  of  using  the  approximate  model  is  very 
much  lower  than  that  of  the  CFD.  (25  seconds  per  function  evaluation,  versus  9000 
seconds.)  Although  the  Taylor  series  (with  10%  move  limits)  improved  upon  the  ini¬ 
tial  design,  it  did  not  find  the  same  design  as  the  CFD,  and  the  increase  m  net  thrust 
was  modest,  at  70%.  Although  the  Taylor  scries  optimization  did  not  exhibit  the 
erratic  behavior  of  the  Gl.A,  it  also  suffers  (to  a  lesser  degree)  the  inability  to  account 
for  the  movement  of  shock  impingement  points.  (Sec  Section  5.5.1.)  The  reason 
that  the  problem  is  less  severe  is  that  the  smearing  of  the  shock,  due  to  tlic  numerical 
dissipation,  conveys  the  information  that  the  shock  is  moving  into  a  new  region,  even 
before  it  gets  there.  In  addition,  the  issue  ol  the  disparity  in  shock-capturing  “crisp¬ 
ness”  between  the  CFD  and  the  approximate  model,  which  hurt  the  GLA,  is  not  an 
issue  in  the  Taylor  .scries  approach,  since  the  approximate  model  is  not  uscc!  there. 
Nevertheless,  the  Taylor  series  optimization  did  not  behave  well. 
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(c)  Details  of  inlet  shock  system  for  opumuni  configuration. 


Rg.  S.IO  (cont’d)  Forebody  configurations  and  shock  structures. 


As  mentioned  in  Section  5.5.2.  a  zeiroth-onk^r  GLA  was  then  applied  to  F„, 
instead  of  to  the  point  properties.  This  technique  yielded  an  optimum  design  very 
close  to  that  found  by  the  CFD.  It  is  important  to  note  that  when  such  a  simple  scal¬ 
ing  is  applied  to  the  integral  measures,  this  is  the  same  as  «impiy  multiplying  all  the 
approximate  method  solutions  (within  one  iteration)  by  a  constant  number.  It  is  not 
surprising,  then,  tiiat  there  should  be  excellent  agreement,  at  least  for  an  uncon¬ 
strained  problem,  when  the  approximate  model  alone  does  such  a  good  Job  This  also 
allowed  the  zeroth  order  integral  measure  GLA  to  use  very  wide  move  limits  (20%), 
wiiich  is  a  major  reason  that  only  1  iteration  was  n^quired  by  this  technique.  This  m.  y 
seem  to  belittle  the  significance  of  applying  a  zeroth-order  scaling  to  the  integral  mea¬ 
sures.  However,  if  it  there  arc  constraints  on  any  integral  measures  in  the  optimiza¬ 
tion,  such  scaling  of  the  constraint  values  would  be  necessary  for  proper  activation 
of  the  constraints.  Therefore,  there  is  value  in  using  such  an  approach  to  the  GLA. 
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when  the  approximate  model  captures  the  general  behavior  of  the  CFD  model  with 
high  fidelity.  The  simplicity  of  implementing  this  technique  is  also  attractive. 

Optimization  histories  of  and  are  shown  in  Fig.  5.11a-c.  The 

horizontal  axes  there  represent  function  evaluations,  meaning  the  number  of  CFD 
calls  for  the  CFD  and  Taylor  series  methods,  and  the  number  of  approximate  model 
calls  for  the  approximate  model  method.  (Note  that  the  values  of  shown  are  as 
calculated  by  the  CFD  solver,  except  for  those  for  the  approximate  model.  'Fhis 
accounts  for  the  apparently  much  higher  F,,.  values  for  this  method.)  Fig.  5. lib 
appears  to  show  that  the  Taylor  series  method  w  as  better  than  the  approximate  model 
for  finding  However,  as  can  be  seen  from  Fig.  5.11c  and  TABLE  5.2,  the 

approximate  method  really  docs  a  better  job  of  replicating  the  behavior  of  the  CFD 
solver. 

It  is  clear  from  these  results  that  the  GLA  and  Taylor  series  method  cannot  be 
applied  to  correcting  point  properties  in  the  forebody/inlet  region  vdth  the  same  ease 
as  was  done  in  the  nozzle  region.  Thus,  the  application  of  GLA’s  to  regions  where 
shocks  dominate  the  flow  requires  a  much  different  treatment  than  that  used  in  current 
structural  optimization  applications. 

5.6  Summary 

After  careful  selection  of  design  variables,  a  successful  optimization  of  the  fore¬ 
body/inlet  region  was  performed  using  the  CFD  alone.  The  approximate  method 
(uncoirected)  was  found  to  optimize  to  a  very  similar  design.  The  Taylor  scries  and 
GLA  (applied  to  point  properties)  did  not  perforra  well,  due  to  their  inability  to 
account  for  the  movement  of  shock  impingement  points.  A  method  of  partially  allevi¬ 
ating  this  problem  has  been  proposed  and  tested  in  one  iteration  of  a  GLA  optimiza¬ 
tion.  This  method  is  based  upon  using  non-diraensiona!ized,  pre-  and  pe.st-shock 
coordinate  systems  for  gradient  and  correction  factor  calculations.  This  method 
appears  to  have  promise,  but  more  work  is  needed  in  this  area.  Additional  optimiza¬ 
tions  were  performed  using  the  GLA  applied  to  integral  measures  (thiust,  drag,  and 

A  zeroth-order  application  of  this  method  was  jjucccssful,  and  holds  some  prom¬ 
ise  for  constrained  optimization  problems  in  which  the  approximate  model  replicates 
well  tlic  general  behavior  of  the  CFD  solver. 


Function  Evaluations 

(a)  Optimization  histories  of  (Note:  valuers  shown  are 

from  CFD  solver,  except  appioximate.  method  values.) 


Fig.  5. 1 1  Optimization  histories  for  forebody/inlet. 
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Function  Evaluations 


(c)  Optimization  histories  of 

Fig.  5.11  (cont’d)  Optimization  histories  for  forebody/inlet. 


Chapter  6 


Conclusions 

6. 1  Constructing  Approximation  Concepts  With 
CFD 

It  was  found  that  approximation  concepts,  as  they  are  currently  employed  in 
structural  optimization,  cannot  be  employed  in  aerodynamic  optimization  problems 
without  modification. 

When  forming  an  aerodynamic  approximation  concept,  special  attention  must  be 
paid  to  the  fact  that  the  number  of  CFD  grid  points  may  change  as  the  grid  clustering 
adapts  to  changes  in  the  design  variables.  This  means  thc.t  the  GI.A  or  Taylor  series 
correction  factors  cannot  be  tied  to  CFD  grid  points.  This  problem  had  previously 
been  identified  (Ref.  3)  but  had  not  heretofore  been  solved. 

This  research  effort  found  a  general  solution  to  this  problem,  by  forming  a  new 
“correction  point”  grid  system,  which  is  independent  of  the  size  and  density  of  the 
CFD  grid.  The  correction  point  grid  system  consists  of  a  very  dense  non-dimension- 
alized  1-D  grid  for  each  wall  surface  of  interest,  as  shown  in  Fig.  2.2.  (A  similar  phi¬ 
losophy  was  applied  to  the  inlet  plane,  for  correction  of  the  flow  properties  there.) 
Whenever  the  CFD  solver  is  run,  the  wall  pressures  are  interpolated  into  the  non- 
dimensionalized  grid  for  each  surface.  Similarly,  the  approximate  model  results  are 
either  interpolated  into  this  grid  system  or,  better  yet,  are  calculated  there  originally, 
when  this  is  practical.  (A  case  where  this  is  not  practical  is  for  the  Method  of  Charac¬ 
teristics,  where  the  MOC  grid  is  set  by  the  MOC  solution,  and  cannot  easily  be  deter¬ 
mined  beforehand.) 

One  special  consideration  inherent  in  constructing  a  GLA  or  Taylor  series  model 
with  CFD  as  the  detailed  model  is  that  the  crispness  of  shock  capturing  may  be  differ¬ 
ent  for  the  CFD  solver  and  for  the  approximate  method.  This  can  hmd  to  erroneous 
correction  factors  in  the  GLA. 

In  addition,  there  are  special  problems  involved  in  using  GLA’s  in  fluid  flow 
regimes  in  which  moving  discontinuous  phenomena  occur  (such  as  shock  waves, 
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transition  to  turbulence,  or  separation).  Such  moving  discontinuities,  if  not  properly 
treated,  can  cause  erratic  behavior  of  the  objective  function. 

Another  special  consideration  when  using  CFD  in  optimizations,  especially  when 
the  gradient  calculations  arc  performed  by  “brute  force”  finite  differences,  is  that  the 
quality  of  the  CFD  solution  must  be  very  high.  Artificial  pressure  spikes  or  excessive 
numerical  dissipation  may  not  only  cause  errors  in  the  values  of  the  objective  func- 
tion,  but  may  also  result  in  the  optimizer  calculating  poor  search  directions,  with  a 
resultant  degradation  in  the  optimization  results.  This  underscores  the  desirability  of 
usmg  analytical  or  quasi-analytical  means  of  finding  the  sensitivity  coefficients. 

It  is  important  that  the  “noise  level”  of  the  CFD  calculations  be  established  by 
examining  the  convergence  histories  of  typical  solutions.  For  meaningful  optimiza¬ 
tions,  the  differences  between  successive  CFD  solutions  must  be  above  the  noise 
level.  This  is  most  critical  during  the  gradient  calculations,  since  small  cnors  in  the 
CFD  solution  there  may  result  in  very  poor  search  directions  being  used. 


6  ,2  Choices  of  Design  Variables 

It  has  long  been  recognized  that  the  selection  of  design  variables  can  have  a  dra¬ 
matic  effect  on  the  efficiency  of  an  optimization  scheme.  In  this  research  it  was 
found,  for  the  forebody/inlet,  that  the  choice  of  design  variables  was  crucial  in  deter¬ 
mining  whether  the  optimizer  selected  a  reasonable  optimum  design. 

Specifically,  the  most  inuitive  design  variable  set  (the  angles  and  lengths  which 
control  the  shock  strengths  and  locations)  did  not  produce  good  results,  because  it  did 
not  allow  for  the  simultaneous  changes  in  some  design  parameters  which  good  design 
practice  dictates.  Several  design  variable  sets  using  only  lengths  were  tried,  and  a  sat¬ 
isfactory  one  was  found.  As  yet,  there  are  no  general  rules  for  determining  the  best 
design  variables. 

The  sensitivity  of  this  optimization  problem  to  the  design  variables  is,  in  large 
part,  due  to  the  fact  hat  arbitrary  shock  locations  were  allowed.  When  only  perfectly- 
placed  shocks  are  cansidered,  there  arc  only  two  design  parameters  (the  nose  and 
ramp  angles),  and  the  optimization  proceeds  smoothly  with  these  selected  as  the 
design  variables. 
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The  nozzleVafterbody  optimizations  (which  did  not  have  to  contend  with  shock 
impingements)  were  found  to  be  insensitive  to  which  of  the  two  candidate  sets  was 
used. 


6.3  Nozzle/Afterbody  (3ptimizations 

Taylor  series  apnroximations  and  two  GLA’s  (one  based  on  the  2-D,  steady 
MOC,  and  one  based  on  isentropic  1-D  flow  theory)  were  successfully  applied  to  the 
nozzle/afterbody  region.  All  of  these  methods  gave  optimum  designs  which  were 
very  close  to  that  found  by  using  the  CFD  alone,  with  a  68%  reduction  in  computa¬ 
tional  cost.  The  successes  in  the  nozzle/afterbody  region  highlight  the  great  potential 
of  approximation  concepts  for  aerodynamic  design  optimization. 

6.4  Forebody/Inlet  OpOmizations 

Optimizations  were  performed  on  the  inlet/fore  body  of  the  2-D  scramjei  vehicle, 
in  which  arbitrary  shock  configurations  were  allowed.  The  objective  function  was  the 
net  thrust.  ITje  optimization  using  CFD  alone  proceeded  without  difficulty.  The 
uncorrected  approximate  method  pertoiraed  well,  optimizing  to  a  design  very  close  to 
that  found  using  the  CFD.  The  GLA  (applied  to  point  properties),  however,  was 
unsuccessful,  owing  to  an  inability  to  aiccount  for  the  motion  of  shocks  as  the  design 
variables  were  changed.  Tliis  resulted  in  poor  search  direction  seiection,  and  in 
erratic  behavior  of  the  objective  function  during  the  optimization  iterations.  The  Tay¬ 
lor  scries  method  did  not  encounter  the  same  erratic  behavior,  but  did  not  optimize  to 
the  same  design  that  the  CFD  and  approximate  methods  did. 

To  alleviate  this  problem,  the  GLA  method  was  then  applied  to  the  integral  mea¬ 
sure  (net  thrust).  Although  such  a  first-order  GLA  was  not  successful,  a  zeroth-order 
GLA  was  found  to  successfully  optimize  to  the  same  design  as  the  CFD. 

To  guide  future  research,  the  detailed  behavior  of  the  point  property  GLA  was 
observed,  and  a  partial  solution  to  the  problems  encountered  by  the  GLA  and  Taylor 
series  was  investigated.  The  basis  of  this  is  to  replace  the  non-dimensionaiized  wall 
coordinate  systems  mentioned  in  Section  6.1  with  non- diraensionalized  pre-  and  post 
shock  coordinate  systems,  as  explained  in  Section  5.5.1  T5ie  gradients  and  correc¬ 
tion  factors  are  calculated  in  these  systeins,  and  the  correction  factors  arc  applied 
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there  before  the  pressures  are  converted  back  into  dimensioned  coordinates,  once  the 
location  of  the  shock  impingement  points  are  known.  This  method  was  found  to 
remove  the  erratic  behavior  of  the  GLA.  However,  there  were  still  problems  with  the 
GLA,  in  that  the  values  of  thrust  and  drag  were  too  high,  and  the  optimizer  treated  the 
new  GLA  almost  exactly  as  it  did  the  original  one.  This  may  have  been  due  in  part  to 
the  fact  that  this  type  of  coordinate  system  was  not  also  applied  to  the  inlet  plane  flow 
conditions,  and  to  the  fact  that  the  CFD  and  approximate  method  do  not  capture 
shocks  equally  crisply,  resulting  in  erroneous  correction  factors  in  the  vicinity  of 
shocks. 

This  concept  of  using  non-dimensionalized  pre-  and  post-discontinuity  coordi¬ 
nate  systems  appears  to  be  a  general  one  for  the  application  of  GLA’s  to  any  problem 
involving  discontinuous  phenomena  which  shift  positions  as  the  design  variables  a:'e 
changed.  Examples  of  such  phenomena  are  shocks,  transition  to  turbulence,  and  sep¬ 
aration. 


6.5  Recommended  Further  Research 

For  the  GLA  approach  to  be  used  in  flow  fields  with  shocks  or  other  shifting  dis¬ 
continuities,  a  method  of  properly  accounting  for  the  motions  of  these  phenomena  in 
tlie  gradient  and  correction  factor  calculations  is  essential.  Research  in  this  area 
should  alv«o  address  the  issue  of  the  difference  in  “crispness”  with  which  the  detailed 
and  approximate  methods  capture  such  discontinuities.  If  such  techniques  can  be 
found,  the  GLA  approach  may  be  made  generally  applicable  in  aerodynamic  optimi¬ 
zation  problems. 

Research  should  be  performed  into  finding  rules  for  determining  which  design 
variable  sets  are  likely  to  provide  the  best  results. 

The  optimizations  performed  in  this  work  had  no  active  constraints  with  regard  to 
the  flow  properties.  Only  geometric  constraints  were  active.  The  behavior  of  the 
GLA  in  the  presence  of  active  “behavior”  constraints  must  be  examined. 

Once  these  issues  have  been  resolved,  attempts  should  be  made  to  apply  the 
GLA  concept  to  more  complicated  problems,  such  as  ootiraizing  a  2-D  scramjet  vchi 
clc  with  chemically-reacting  internal  and  external  flows.  3-D  vehicle  optimizations 
present  a  more  challenging  test  of  the  method,  but  also  one  in  which  the  potential  pay¬ 
off  of  the  GLA  is  great. 
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Appendix  A 


The  Metliod  of  Feasible  Directions 


A.  1  Introduction 

Each  scramjet  vehicle  optimization  was  pcrfonned  by  combining  the  well-proven 
optimization  software  package  “Design  Optimization  Tools”,  or  DOT  (See  Ref,  12.), 
with  one  or  more  2-D  flow  solvers  (or  “analysis  programs”).  DOT  uses  the  “Modified 
Method  of  Feasible  Directions”  to  solve  constrained  nonlinear  optimization  problems. 
In  this  case,  DOT  was  instructed  to  find  the  vehicle  design  which  maximized  net 
thrust  (Thrust  -  Drag),  and  which  also  met  certain  constraints  on  vehicle  dimensions 
and  flow  properties  at  critical  points  on  the  vehicle.  The  flow  solvers  provided  DOT 
with  the  values  of  net  thrust  and  of  the  flow  properties  considered  in  the  constraints. 
This  chapter  describes  the  salient  features  of  the  Method  of  Feasible  Directions. 

A.  2  Optimization  Method 


The  general  optimization  problem  can  be  formulated  as  follows.  Let  the  objective 
function  (in  this  problem,  the  net  thrust)  be  F(A0,  where  X  is  the  vector  of  design  vari¬ 
ables  (e.g.  nose  angle,  inlet  compression  ramp  angle,  etc.).  The  elements  of  X  are  X] 
X2 ,  etc.  Let  g  be  the  vector  of  inequality  constraints  (e.g.  pressure  at  the  end  of  the 
nozzle  wall  must  be  greater  than  some  value,  to  prevent  reverse  flow  in  the  nozzle). 
Define  h  as  the  vector  of  equality  constraints.  Furthermore,  let  the  elements  of  X  be 
limited  to  certain  ranges  of  values  (side  constraints).  Following  the  terminology  in 
Ref.  11,  we  now  write  the  problem  statement  as: 


Maximize  F(X) 
Subject  to: 

/»*(X)  =0 
XjSX.sXr 


Objective  function 
j  =  l,m  Inequality  constraints 
k  =  1,/  Equality  constraints 
i  =  l,n  Side  constraints 
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where 


X  = 


fX. 


Design  variables 


The  vector  X  defines  a  design  space.  Each  point  in  the  design  space  has  associ¬ 
ated  with  it  a  value  of  the  objective  function.  The  constraints  can  be  thought  of  as 
boundaries  within  the  space.  Fig.  A.l  portrays  such  a  problem  which  has  only  two 
design  variables.  (This  figure  is  based  on  Fig.  1.6  of  Ref.  11.)  The  basic  concept  of 
the  Method  of  Feasible  Directions  is  that,  as  we  search  the  design  space  for  an  opti¬ 
mum,  we  search  in  directions  which  both  improve  the  value  of  the  objective  function, 
and  do  not  violate  any  constraints.  At  point  A  in  Fig.  A.l  the  constraint  is  active.  If 
we  try  to  search  in  the  direction  of  VgiX) ,  we  will  only  increase  the  violation  of  the 
constraint.  The  only  way  to  guarantee  that  the  constraint  will  not  be  violated  (at  least 
for  a  small  step)  is  to  limit  the  search  to  directions  which  lie  on  the  opposite  side  of  a 
Unc  normal  to  Vg  (X)  ,  All  such  directions  are  defined  as  being  “feasible”  directions. 
Similarly,  the  only  way  to  guarantee,  for  a  small  step,  that  moving  in  a  given  search 
direction  will  result  in  an  increase  in  the  objective  function,  is  to  require  that  direction 
to  lie  in  the  half-plane  which  is  bounded  by  the  normal  to  VF  (X) ,  and  which  includes 
the  direction  VF  (X) ,  All  such  directions  are  “usable”  directions.  In  snort,  we  should 
search  in  a  direction  which  is  to  some  degree  “along”  V  F  (X)  and  “opposite  to” 
Vg  (X) .  Or,  we  seek  a  direction  which  is  both  “usable”  and  “feasible”.  The  usable- 
feasible  sector  at  point  A  in  Fig,  A.  1  contains  all  the  usable-feasible  directions  ema¬ 
nating  from  that  point.  Note,  however,  that  as  we  move  in  such  a  directicn,  we  may 
soon  find  ourselves  in  violation  of  the  constraint.  The  usability  and  feasibility  of  a 
direction  is  guaranteed  only  in  the  immediate  vicinity  of  the  originating  point. 

At  point  B,  there  is  no  usable-feasible  direction,  so  this  point  is  the  optimum 
design  for  this  simple  example.  Fhe  optimality  criterion  in  this  case  is  that  V  F  (X) 
and  Vg  (X)  are  in  exactly  opposite  directions. 
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Fig.  A.1  Example  of  two- variable  design  space  with 
inequality  constraint. 


To  determine  the  optimality  criteria  for  a  more  general  design  problem,  we  again 
consider  the  gradients  of  the  objective  function  and  constraints,  as  shown  in  Fig.  A.2a. 
Point  X*  is  the  constrained  optimal  point  in  this  space.  That  point  is  nestled  in  the 
“valley”  formed  by  the  inequality  constraints.  This  can  be  expressed  mathematically 
by  noting  that,  for  such  a  configuration,  the  gradients  of  the  two  constraints  straddle 
VFiX) .  Thus  we  could  construct  a  linear  combination  of  Vg,  (^X)  and  V  g2  (JQ  which 
will  be  exactly  equal  to  V  F  (X) .  This  way  of  denoting  that  the  optimal  point  is  at  the 


(a)  Objective  function  and  constraint  geometry. 


Feasible  Sector 

wit  ^2 


(b)  Usable  and  Feasible  Directions  at  optimal  point. 


Fig.  A.2  Illustration  of  Kuhn-Tucker  optimality  conditions 
for  a  two-design  variable  problem  with  two  active  constraints. 
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bottom  of  the  valley  is  embodied  in  the  well-known  Kuhn-Tucker  necessary  condi¬ 
tions  for  optimality. 

1.  X*  is  feasible.  (A.l) 

2.  Xjgj(X)=o  X.5  0  (A.2) 

m  t 

3.  +  +  =0  (A.3) 

jml  i.l 

Xj^O 

unrestricted  in  sign. 

Equations  (A.l)  and  (A.2)  simply  require  that  the  design  be  feasible  and  satisfy  all 
inequality  constraints.  Eq.  (A.3)  requires  that  some  linear  combination  of  the  gradi¬ 
ents  of  the  active  constraints  be  equal  to  V  F  (AO .  As  shown  in  Fig.  A.2b,  this 
implies  that  there  is  no  usable-feasible  search  direction  at  this  point.  Note  that  both 
inequality  and  equality  constraints  are  considered  in  Eq.  (A.3).  Note  also  that  if  no 
constraints  are  active,  then  Eq.  (A.3)  says  that  the  optimal  design  occurs  where 
V  F(X)  -  0.  This  is  analogous  to  the  familiar  “zero  slope”  condition  for  the  extre¬ 
mum  (maximum,  minimum  or  inflection  point)  of  a  function  of  a  single  variable. 

The  Kuhn-Tucker  conditions  can  only  guarantee  that  a  local  extremum  has  been 
reached.  They  could  just  as  easily  indicate  a  local  minimum  or  saddle  point  as  they 
would  a  local  maximum.  The  search  technique  used  will  help  avoid  settling  on  the 
wrong  type  of  extremum.  A  further  difficulty  is  that  there  may  well  be  other  local 
optima  which  yield  higher  values  of  F(X).  For  problems  which  are  highly  non-linear, 
it  is  usually  necessary  to  choose  an  initial  design  in  the  neighborhood  of  the  global 
optimum  to  guarantee  that  the  global  optimum  will  be  reached.  However,  even  if  the 
global  best  is  not  found,  this  method  will  usually  produce  a  design  which  is  better 
than  the  initial  design.  For  many  engineering  applications,  this  may  be  sufficient  jus¬ 
tification  for  using  such  techniques. 

The  first  step  in  implementing  the  Method  of  Feasible  Directions  with  DOT  is  to 
select  an  initial  design,  X°.  Then  DOT  calculates  an  initial  search  direction,  which  is 
in  the  direction  of  V  F(X),  unless  the  initial  design  is  infeasible.  In  that  case,  DOT 
first  seeks  to  relieve  the  constraint  violation,  while  improving  the  value  of  F(X),  if  this 
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is  possible.  All  evaluation.*;  of  F(X),  V  F(X)  ,  and  V  g(X)  are  provided  by  a  sepa¬ 
rate  “analysis  program”  (in  this  research,  the  flow  solver  program).  If  the  analysis 
program  does  not  provide  the  gradients  directly,  then  they  are  calculated  by  repeated 
evaluations  of  F(X)  and  g(X),  wherein  the  elements  of  X  are  perturbed  about  the  ini¬ 
tial  design,  X°.  This  finite-difference  method  of  finding  the  elements  of  VF(;ir) 
(which  are  commonly  referred  to  as  the  “sensitivity  coefficients”  of  F(X)  )  can  be 
extremely  costly  in  computer  resources.  Much  work  has  been  done  in  recent  years  on 
determining  aerodynamic  sensitivity  coefficients  analytically  or  by  approximate 
numerical  techniques,  to  avoid  having  to  use  finite  differencing.  (See  Ref’s.  4,  13,  and 
14,  for  example.)  However,  the  present  work  used  finite  differences  exclusively. 

A  “one-dimensional  search”  is  performed  along  the  search  direction,  until  the 
optimal  constrained  value  of  V  F  (X)  for  this  direction  is  found.  Then  a  new  search 
direction  is  calculated.  This  direction  is  not  in  the  direction  of  V  F(X) ,  but  is  rather 
in  a  direction  determined  by  the  Conjugate  Direction  Method  or  the  Variable  Metric 
Method.  (See  Ref.  11.)  Such  methods  achieve  much  faster  convergence  than  using 
the  gradient  direction.  As  the  search  approaches  a  constraint,  “push-off  factors”  are 
used  in  the  constraint  calculations  to  force  the  search  back  into  the  feasible  region, 
before  the  constraint  is  actually  violated.  A  new  one-dimensional  search  is  per¬ 
formed,  and  the  process  continues  until  no  further  improvement  can  be  made  in  the 
objective  function  without  violating  any  constraint.  The  resulting  design  is  then  con¬ 
sidered  to  be  optimal. 

The  major  limitations  of  this  type  of  optimization  scheme  are:  that  it  requires  that 
the  objective  function  and  constraint  functions  be  continuous;  and  that  it  will  only 
guarantee  a  local  optimum  will  be  found.  Since  many  engineering  problems  do  not 
have  continuous  objective  and  constraint  functions,  the  first  limitation  would  appear 
to  severely  limit  the  usefulness  of  the  method.  But  is  has  been  found  in  practice  that 
some  degree  of  discontinuity  can  be  tolerated.  The  problem  of  determining  if  the 
optimal  solution  is  actually  a  global  optimum  can  be  somewhat  alleviated  by  running 
the  optimization  several  times  starting  at  different  initial  designs.  This  may  be  quite 
expensive,  particularly  if  there  is  a  large  number  of  design  variables. 


91 


A.  3  Summary 

The  optimization  technique  used  in  this  research  employs  the  Modified  Method 
of  Feasible  Directions,  as  embodied  in  the  optimization  package  DOT.  This  scheme 
requires  function  and  constraint  evaluations  and  their  gradients  (with  respect  to  the 
design  variables),  which  arc  provided  1  .  or  nriore  analysis  programs.  The  gradi¬ 

ents  were  found  using  the  “bnite  force’  jnite-difference  method.  The  optimization 
method  establishes  a  search  direction  .i  the  usable-feasible  sector  (if  one  exists)  and 
then  performs  a  1 -dimensional  sea.,,11  in  this  direction.  When  the  optimal  point  in 
this  direction  is  found,  a  new  search  direction  is  determined  and  a  new  1 -dimensional 
search  is  conducted.  The  process  is  repeated  until  convergence  is  reached.  The 
Kuhn-Tucker  necessary  conditions  for  optimality  are  used.  These  conditions  can  only 
guarantee  that  a  local  optimum  is  reached. 
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